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AEROELASTIC, CFD AND DYNAMICS COMPUTATION AND 
OPTIMIZATION FOR BUFFET AND FLUTTER APPLICATIONS 

GRANT NO. NAG-1-648 
Osama A. Kandil 

Aerospace Engineering Department 
Old Dominion University, Norfolk, VA 23529-0247 


Accomplishments 

In the period of December 1, 1996 to November 30, 1997, the Principal Investigator (PI) 

along with the assistance of two Ph.D. students have achieved the following accomplishments 

under this grant: 

I. Publications 

1.1 Kandil, O. A. and Sheta, E. F., “Coupled and Uncoupled Bending-Torsion Responses of 

Twin-Tail Buffet,” Accepted for publication in Journal of Fluids and Structures, Academic 
Press, to appear 1998. ^ 

1.2 Kandil, O. A. and Sheta, E. F., “Coupled and Uncoupled Bending Torsion Responses of 
Twin Tail Buffet,” Fourth International Symposium on Fluid/Structure, Aeroelasticity and 
Flow Induced Vibrations and Noise, ASME 97-129. 1997 ASME International Engineering 
Congress and Exposition , Dallas, Texas, November 1997, AD-Vol . 53-3, pp. 1-12 (A 
copy is attached) . 

1.3 Kandil, 0. A., Liu, C. H. and Sheta, E. F., “Effects of Coupled and Uncoupled Mode 
Responses on Tail Buffet Over a Wide Range of Angles of Attack,” IUTAM- Symposium on 
Slender Vortex Dynamic, Aerodynamics Institute, Aachen, Germany , August 31 -September 
4, 1997. Invited keynote paper . 

1.4 Kandil, 0. A., and Abdelhamid, Y. A., “Computation and Validation of Delta Wing 
Pitching up to 90° Amplitude,” A LA A 97-3573-CP, AIAA Atmospheri c Flight Mechanics 
Conference, New Orleans, LA, August 11-13, 1997, pp. 221-231. 

1.5 Kandil, 0. A., Sheta, E. F., Massey, S. J., “Fluid/Structure Twin Tail Buffet Responses 
Over a Wide Range of Angles of Attack,” AIAA 97—2261, AIAA Applied Ae rodynamics 
Conference, Atlanta, GA , June 23-25, 1997, pp. 373—386. (A copy i s attached). 

1.6 Kandil, 0. A. and Menzies, M. A., “Effective Control of computationally Simulated Wing 
Rock in Subsonic Flow,” AIAA 97-0831 , AIAA 35th ASM, Reno, Nevada, January 1997. 

1.7 Kandil, 0. A., “Recent Advances in Multidisciplinary Aeronautical Problems of 
Fluid/Structures/Dynamics Interaction,” International Seminar Series II Proceedings^ 
Institute of Aeronautics and Applied Mechanics, Warsaw Univers ity of Technology, 
November 1996. Invited paper , pp. 15-34. (A copy is attached) . 




II. Papers Accepted for Conference Presentation 

II I- E ssam ’ E. P* Kandil, O. A., “Effect of Configuration Pitching Motion on Twin Tail 
Buffet Response,” AIAA 98-0520 . AIAA 36th ASM, Reno, NV, January 1998. 

11.2. Abdelhamid, Y. A. and Kandil, O. A., “Effect of Reduced Frequency on Super-Manuever 
Delta Wing,” AIAA 98-0414 . AIAA 36th ASM, Reno, NV, January 1998. 

11.3. Massey, S. J. and Kandil, 0. A., “Effect of Apex Flap Deflection on Vertical Tail Buffet,” 
AIAA 98—0762 , AIAA 36th ASM Reno, NV, January 1998. 


III. Proposal Submitted for CRAY-C-90 Usage 

A proposal has been submitted for usage of the National Aerodynamic Simulation (NAS) 
Facilities at NASA Ames Research Center, CA, on May 1, 1997. A total of 200 hrs. of C-90 
CPU time has been requested. 

A technical summary has also been submitted to NAS on November 26, 1997. It summarized 
the achievements, goals, approach and recent results of effects of coupled and uncoupled bending 
and torsion buffet responses of a typical twin-tail configuration. Comparison of these results is 
shown with the available experimental data. __ 


IV. Graduate Students 

Two Ph.D. students have been assisting the P,I. to carry out the tasks of this grant and write 
their Ph.D. dissertations. The status of the students is given below: 


VI. Mr. Steven J. Massey (U.S. Citizen): He has been working on his Ph.D. degree since 
May 1994, and has been supported under this grant and a fellowship from the Aerospace 
Engineering Department. He successfully defended his Ph.D. dissertation in October, 1997 
and is graduating this semester. His Ph.D. dissertation is titled, “Development of A Coupled 
Fluid/Structure Aeroelastic Solver with Applications to Vortex Breakdown Induced Twin Tail 
Buffet.” Drs. Woodrow Whitlow, Jr. and Robert Bennett served on his Ph.D. committee. A 
copy of his dissertation abstract is attached . 

VII. Mr. Mark W. Flanagan (U.S. Citizen): He has been working on his Ph.D. degree' since 
January 1994. He has been supported under a fellowship from the AE Dept, and work at 
Dynamic Engineering Inc. Currently, he is a Ph.D. candidate (since Spring 1996) and he 
is expected to finish his Ph.D. degree in December 1998. His Ph.D. dissertation focuses on 
simulation and optimization control of tail buffet in supersonic internal vortex breakdown 
flows in a configured duct. Drs. Woodrow Whitlow, Jr. and Robert Bennett are serving 
on his Ph.D. committee. 








A third Student, Mr, Essam Sheta, who is supported by the Aerospace Engineering 

Department, is also writing his Ph.D. Dissertation on Flow Turbulence and Control Effects 

on Coupled and Uncoupled Bending Torsion Modes of Twin Tail Response. He is expected to 

finish his Ph.D. dissertation in May 1998. 

V. Technical Group Participations and Conference Presentations: The past year 

(12/1/96-11/30/97) was a very productive one for technical group participations and 
conference presentations by the P.I. and his students. These are listed below: 

1. “Coupled and Uncoupled Bending torsion Responses of Twin Tail Buffet," ASME 97-129, 
1997 ASME International Engineering Congress and Exposition, Dallas, Texas , November 
1997. 

2. “Effects of Coupled and Uncoupled Modes on Tail Buffet Over a Wide Range of Angles of 
Attack,” rUTAM Colloquim on Slender Vortex Dynamics, Aerodynamics Institute, A achen. 
Germany , August 31 — September 3, 1997. (Invited) 

3. “Computation and Validation of Delta wing Pitching up to 90° Amplitude,” ALAA Atmo- 
spheric Flight Mechanics Conference, New Orleans, LA, August 11-13, 1997. 

4. “Fluid/Structure Twin Tail Buffet Responses Over a Wide Range ofAngles of Attack,” AIAA 
Applied Aerodynamics Conference, Atlanta, GA, June 23-25, 1997. 

5. “Effect of Turbulence on Twin Tail-Buffet Response,” VCES Spring Conference , Hampton, 
VA , June 10, 1997. 

6. “Computation and Validation of Delta Wing Pitching Response for very high Angles of 
Attack,” VCES Spring Conference, Hampton, VA, June 10, 1997. 

7. “Effective Control of computationally Simulated Wing Rock in Subsonic Flow,” AIAA 35th 
Aerospace Sciences Meeting, Reno, Nevada, January 9—12, 1997. 

8. General Chairman and Organizer of AIAA Applied Aerodynamics (APA) Program. 
Aerospace Sciences Meeting, Reno, NV, January 1998. Two new special sessions which 
are related to the current grant work have been introduced in the APA program for the first 
time. These are Aero/Structures Interaction and Flow control using oscillating jets. 

9. The P.I. was invited by the Ohio Aerospace Institute to serve on the National Committee 
of Fixed Wng Vehicle Initiative for Aerodynamics, Structures and Control Integration 
(ASCI), initiated by Industry/Govemment/Academia Technical Leaders, also a member of 
the Unsteady Flow Team and Academia Team, May-July 1997. A major report has been 
produced for the next 15 year research needs in this important area. 

10. The P.I. is planning to participate in the coming workshop at NAS A Langley Research on the 
Nonlinear Aeroservoelasticity, organized by the Branch in the period of January 20-22, 1998. 
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ABSTRACT * 

The effect of coupled and uncoupled bend- 
ing and torsion modes on flexible twin-tail buf- 
fet is considered. This multidisciplinary problem 
is investigated using three sets of equations on a 
multi-block grid structure. The first set is the un- 
steady, compressible, full Navier-Stokes equations 
which are used for obtaining the flow-filed vec- 
tor and the aerodynamic loads on the twin tails. 
The second set is the coupled aeroelastic equations 
which are used for obtaining the bending and tor- 
sional deflections of the twin tails. The third set 
is the grid-displacement equations which are used 
for updating the grid coordinates due to the tail 
deflections. The configuration is pitched at 30° 
angle of attack and the freestream Mach number 
and Reynolds number are 0.3 and 1.25 million, re- 
spectively. Keeping the twin tails as rigid surfaces, 
the problem is solved for the initial flow condi- 
tions. Next, the problem is solved for the flexible 
twin tails responses due to the unsteady loads pro- 
duced by the vortex breakdown flow of the delta- 
wing leading-edge vortex cores. The configura- 
tion is investigated for the effect of coupled and 
uncoupled bending and torsion modes using two 
different separation distances of the twin-tail; the 
inboard and the outboard positions. 

1 INTRODUCTION 

In order to maximize the effectiveness of the 
fighter aircraft that operate well beyond the buffet 
onset boundary, the design of the new generation 
of fighter aircraft should account for both high 
maneuver capabilities and the aeroelastic buffet 
characteristics at high and wide range of angles of 
attack. The maneuver capabilities are achieved, 
for example in the F/A-18 fighter, through the 


combination of the leading-edge extension (LEX) 
with a delta wing and the use of vertical tails. 
The LEX maintains lift at high angles of attack 
by generating a pair of vortices that trail aft over 
the top of the aircraft. The vortex entrains air 
over the vertical tails to maintain stability of the 
aircraft. At some flight conditions, the vortices 
emanating from the highly-swept LEX of the delta 
wing breakdown before reaching the vertical tails 
which get bathed in a wake of unsteady highly- 
turbulent, swirling flow. The vortex-breakdown 
flow produces unsteady, unbalanced loads on the 
vertical tails and causes a peak in the pressure 
spectrum that may be tuned to different struc- 
tural modes depending on the angle of attack and 
dynamic pressure. This in turn produces severe 
buffet on the tails and has led to their premature 
fatigue failure. If the power spectrum of the tur- 
bulence is accurately predicted, the intensity of 
the buffeting motion can be computed and the 
structural components of the aircraft can be de- 
signed accordingly. 

Experimental investigation of the vertical 
tail buffet of the F/A-18 models have been con- 
ducted by several investigators such as Sellers, et 
al. (1988), Erickson, et al. (1989), Wentz (1987), 
Lee and Brown (1990), and Cole, et al. (1990). 
These experiments showed that the vortex pro- 
duced by the LEX of the wing breaks down ahead 
of the vertical tails at angles of attack of 25° 
and higher producing unsteady loads on the ver- 
tical tails, and the buffet response occurs in the 
first bending mode, increases with increasing dy- 
namic pressure and is larger at M — 0.3 than that 
at higher Mach numbers. Bean and Lee (1994) 
showed that buffeting in the torsional mode oc- 
curred at a lower angle of attack and at larger lev- 
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els compared to the fundamental bending mode. 
An extensive experimental investigation has been 
conducted to study vortex-tail interaction on a 
76° sharp-edged delta wing with vertical twin-tail 
configuration by Washburn, Jenkins and Ferman 
(1993). The vertical tails were placed at nine 
locations behind the wing. The experimental 
data showed that the aerodynamic loads are more 
sensitive to the chordwise tail location than its 
spanwise location. As the tails were moved later- 
ally toward the vortex core, the buffeting response 
and excitation were reduced. 


steady, compressible, Reynolds- averaged Navier- 
Stokes equations. The second set consists of 
the aeroelastic equations for coupled bending and 
torsional modes. For uncoupled bending-torsion 
modes, x$ is set equals to zero in Eqs (4) and 
(5). The third set consists of equations for de- 
forming the grid according to the twin tail deflec- 
tions. Next, the governing equations of each set 
along with the initial and boundary conditions are 
given. 

2,1 Fluid-Flow Equations: 


Kandil, Kandil and Massey (1993) pre- 
sented the first successful computational simula- 
tion of the vertical tail buffet using a delta wing- 
single flexible vertical tail configuration. The tail 
was allowed to oscillate in bending modes. Un- 
steady vortex breakdown of leading-edge vortex 
cores was captured, and unsteady pressure forces 
were obtained on the tail. Later on, Kandil, et al. 
(1994-1995) allowed the vertical tail to oscillate in 
both bending and torsional modes. The total de- 
flections and frequencies of deflections and loads 
of the coupled bending-torsion case were found to 
be one order of magnitude higher than those of the 
bending case only. Also, it lias been shown that 
the tail oscillations change the vortex breakdown 
locations and the unsteady aerodynamic loads on 
the wing and tail. 

The buffet responses of twin-tail model has 
been studied by Kandil, Sheta and Liu (1996). 
The twin tails were considered at a = 30° and for 
three different spanwise positions of the twin tails. 
A multi-block grid structure was used to solve the 
problem. The loads, deflections, frequencies and 
root bending moments were reduced as the twin 
tails moved laterally toward the vortex core. The 
outboard position of the tails produced the least 
of these responses. In a recent paper by Kandil, 
Sheta and Massey (1997), the buffet response of 
twin-tail model in turbulent flow was considered 
at a wide range of angles of attack. The computa- 
tional results were in good quantative agreement 
with the experimental data of Washburn, et al 
(1993). 

In this paper, we consider the effect of cou- 
pled and uncoupled bending and torsion modes on 
the flexible twdn-tail buffet response for two differ- 
ent spanwise separation distance of the twin- tail; 
the inboard position (33% wing span) and the 
outboard position (78% wing span). 


2 FORMULATION 

The formulation consists of three sets of 
governing equations along with certain initial and 
boundary conditions. The first set is the un- 


The conservative form of the dimensionless, 
unsteady, compressible, full Navier-Stokes equa- 
tions in terms of time-dependent, body-conformed 
coordinates £ 2 and £ 3 is given by 


dQ 0 Em 
dt + d£ m 


d(E v ) s 

d£ s 


= 0; m = 1 -3, s 


1-3 (1) 


where 


C 1 = ^ 7T, (xi,j 2 i ^3,0 (2) 

Q = j[p,pu x ,pu 2 ,pu3,pe]\ (3) 

E m and {E v ) s are the £ m -inviscid flux and 
£ 5 -viscous and heat conduction flux, respectively. 
Details of these fluxes are given by Kandil, Kandil 
and Massey (1993). 


2.2 Aeroelastic Equations: 

The dimensionless, linearized governing equa- 
tions for the coupled bending and torsional vibra- 
tions of a vertical tail that is treated as a can- 
tilevered beam are considered. The tail bending 
and torsional deflections occur about an elastic 
axis that is displaced from the inertial axis. These 
equations for the bending deflection, w r , and the 
twist angle, 6 , are given by 


d 2 

dz 2 


E n-) 

8 2 9 


d 2 w . 


+m(z)x e (z) w (z,t) = N(z,t) 


0 

dz 


GJ(z), 


89 


dz 


- m{z)x e (z)^Z-(z,t) 


dt 2 


d 2 9 




(4) 


(5) 


where z is the vertical distance from the 
fixed support along the tail length, /(, El and GJ 
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the bending and torsional stiffness of the tail sec- 
tion, m the mass per unit length, I $ the mass- 
moment of inertia per unit length about the elas- 
tic axis, N the normal force per unit length, M t 
the twisting moment per unit length and xq the 
distance between the elastic axis and inertia axis. 
When xe = 0.0 the bending and torsion modes 
are dynamically decoupled. The characteristic pa- 
rameters for the dimensionless equations are c*, 
a’ 0 , pic and c*/a^ for the length, speed, den- 
sity and time; where c* is the delta wing root- 
chord length, the freestream speed of sound 
and p ^ the freestream air density. The details 
of the solution method are given in Kandil, Sheta 
and Massey (1997). 


2.3 Grid Displacement Equations: 

Once w and 6 are obtained at the n 4- 1 time 
step, the new grid coordinates are obtained us- 
ing simple interpolation equations. In these equa- 
tions, the twin tail bending displacements, , 
and their displacement through the torsion angle, 
are interpolated through cosine functions. 


2.4 Boundary and Initial Conditions: 

Boundary conditions consist of conditions 
for the fluid flow and conditions for the aeroelas- 
tic bending and torsional deflections of the twin 
tail. For the fluid flow, the Riemann-invariant 
boundary conditions are enforced at the inflow 
and outflow boundaries of the computational do- 
main. At the plane of geometric symmetry, pe- 
riodic boundary conditions are specified. On the 
wing surface, the no-slip and no-penetration con- 
ditions are enforced and = 0. On the tail 
surface, the no-slip and no-penetration conditions 
for the relative velocity components are enforced 
(points on the tail surface are moving). The nor- 
mal pressure gradient is no longer equal to zero 
due to the acceleration of the grid points on the 
tail surface. This equation becomes = -pa t .h, 
where a* is the acceleration of a point on the tail 
and n is the unit normal. 

The initial conditions of the fluid flow cor- 
respond to the freestream conditions with no- 
slip and no-penetration conditions on the wing 
and tail. For the aeroelastic deflections of the 
tail, the initial conditions for any point on the 
tail are that the displacement and velocity are 
zero, w(z, 0) = 0, §f(*,0) = 0, 0(z, 0) = 0 and 

ff( 2 ’0) = 0. 


3 METHOD OF SOLUTION 

The first step is to solve for the fluid flow 
problem using the vortex- breakdown conditions 
and keeping the tail as a rigid beam. Navier- 
Stokes equations are solved using the implicit, 
flux-difference splitting finite-volume scheme. The 
grid speed is set equal to zero in this step. 
This step provides the flow field solution along 
with the pressure differences across the tails. The 
pressure differences are used to generate the nor- 
mal force and twisting moment per unit length 
of each tail. Next, the aeroelastic equations are 
used to obtain the twin tail deflections, ivi,j % k and 
0i 3 k- The grid displacement equations are then 
used to compute the new grid coordinates. The 
metric coefficient of the coordinate Jacobian m$- 
trix are updated as well as the grid speed, -|j-. 
This computational cycle is repeated every time 
step. 

4 COMPUTATIONAL APPLICATIONS 
AND DISCUSSION 

4.1 Twin Tail-Delta Wing Configuration: 

The twin tail-delta wing configuration con- 
sists of a 76°-swept back, sharp-edged delta wing 
(aspect ratio of one) and dynamically scaled flex- 
ible twin tails similar to those used by Washburn, 
et. al. (1993). The vertical tails are oriented 
normal to the upper surface of the delta wing and 
have a centerline sweep of 53.5°. A multi- block 
grid consisting of 4 blocks is used for the solution 
of the problem. The first block is a O-H grid for 
the wing and upstream region, with 101X50X54 
grid points in the wrap around, normal and axial 
directions, respectively. The second block is a H-H 
grid for the inboard region of the twin tails, with 
23X50X13 grid points in the wrap around, nor- 
mal and axial directions, respectively. The third 
block is a H-H grid for the outboard region of the 
twin tails, with 79X50X13 grid points in the wrap 
around, normal and axial directions, respectively. 
The fourth block is a O-H grid for the down- 
stream region of the twin tails, with 101X50X25 
grid points in the wrap around, normal and axial 
directions, respectively. Figure 1 shows the three 
dimensional grid topology and a front view blow- 
up of the twin tail-delta wing configuration. 

Each tail is made of a single Aluminum spar 
and Balsa wood covering. The Aluminum spar has 
a taper ratio of 0.3 and a constant thickness of 
0.001736. The chord length at the root is 0.03889 
and at the tip is 0.011667, with a span length 
of 0.2223. The Aluminum spar is constructed 
from 6061-T6 alloy with density, p. modulii of 
elasticity and rigidity, E and G, of 2693 kg/m 3 , 
6.896X10 10 N/m 2 and 2.5925X10 10 N/m 2 ; re- 
spectively. The corresponding dimensionless 
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quantities are 2198, 4.595 x 10 5 and 1.727 x 10 5 ; 
respectively. 

The Balsa wood covering has a taper ra- 
tio of 0.23 and aspect ratio of 1.4. The chord 
length at the root is 0.2527 and at the tip is 
0.058, with a span length of 0.2223. The Balsa 
thickness decreases gradually from 0.0211 at the 
tail root to 0.0111 at the tail midspan and then 
constant thickness of 0.0111 is maintained to 
the tail tip. The tail cross section is a semi- 
diamond shape with bevel angle of 20°. The 
Balsa density, modulii of elasticity and rigidity, E 
and G, are 179.7 kg/m 3 , 6.896AM0 8 N/m 2 and 
2.5925AM 0 8 N/m 2 ] respectively. The correspond- 
ing dimensionless quantities are 147, 4.595 x 10 3 
and 1.727 x 10 3 ; respectively. The tails are as- 
sumed to be magnetically suspended and the lead- 
ing edge of the tail root is positioned at x/c = 1.0, 
measured from the wing apex. The configura- 
tion is pitched at 30° angle of attack and the 
freestream Mach number and Reynolds number 
are 0.3 and 1.25 x 10 6 ; respectively. The config- 
uration is investigated for two spanwise positions 
of the twin tails; the inboard position and the 
outboard position corresponding to a separation 
distance between the twin tails of 33% and 78% 
of the wing span; respectively. 

Keeping the twin tail as rigid surfaces, the 
unsteady full Navier-Stokes equations are inte- 
grated time accurately using the implicit, flux- 
difference splitting scheme of Roe with Reynolds 
number of 1.25 million and angle of attack of 30°. 
The initial conditions are obtained after 10,000 
time steps with At — 0.001. Next, the results 
of the coupled and uncoupled bending and tor- 
sion modes are presented. For the coupled bend- 
ing and torsion case, the inertia axis is assumed 
downstream the elastic axis by xq = 0.003. For 
the uncoupled case, x$ — 0.0. 

4.2 Uncoupled Bending- Torsion Modes: 

Figure 2 shows three-dimensional views of 
the leading-edge vortex cores particle traces and 
iso-total pressure surfaces. Figures 3 and 4 show 
front views for the total pressure contours on the 
wing surface and in cross flow planes at x = 1.03 
and x = 1.22, and the instantaneous streamlines 
of the inboard tail position. The leading-edge 
vortex cores experience asymmetric breakdown 
upstream of the twin tail due to the upstream 
effect of the tails motion. The vortices are totally 
outboard of the twin tail. The cores are moved 
upward as the flow traveled downstream. Smaller 
size vortex cores appear underneath the primary 
wing vortex and it becomes larger in size as it 
travels downstream. These are the tail vortices 
observed by Washburn (1993). The tail vortices 


exist at the outer surfaces of the tails and they are 
rotating in the opposite direction to those of the 
primary wing vortices. Figure 5 shows the distri- 
bution of the surface pressure coefficient covering 
the wing from x = 0.3 tor = 1.0. The largest suc- 
tion peaks are pronounced at the position of the 
wing vortex cores, and the peaks values decreased 
as the flow travels downstream. Figure 6 shows 
the distribution of the leading-edge total struc- 
tural deflection and the root bending moment for 
the left and right tails for 20 dimensionless time 
after the initial conditions. The tails deflections 
are in first, second and third mode shapes. The 
frequencies of the torsion deflections are almost 
twice those of the bending deflections. 

Figure 7 shows three-dimensional views of 
the leading-edge vortex cores particle traces and 
iso-total pressure surfaces. Figures 8 and 9 show 
front views for the total pressure contours on the 
wing surface and in cross flow planes at x = 1.03 
and x = 1.22, and the instantaneous stream- 
lines of the outboard tail position. The tails cut 
through the vortex breakdown flow of the leading- 
edge vortex cores. The tail vortices are also out- 
board of the tails. The location of the vortex core 
with respect to the tail produces an increase in 
the aerodynamic damping, causing the tail de- 
flection to decrease. The tail vortices are also 
shown to rotate in the opposite direction to those 
of the primary wing vortices. Figure 10 shows 
the distribution of the surface pressure coefficient 
covering the wing from x = 0.3 to x = 1.0. The 
suction peaks are less than those of the inboard 
twin-tail case. Figure 11 shows the distribution 
of the leading-edge total structural deflection and 
the root bending moment for the left and right 
tails for 20 dimensionless time after the initial 
conditions. The levels of loads and defections are 
much lower than those of the inboard twin-tail 
case. The tails are shown to oscillate in one di- 
rection only in first and second mode shapes. 


4.3 Coupled Bending- Torsion Modes: 

Figure 12 shows three-dimensional views of 
the leading-edge vortex cores particle traces and 
iso-total pressure surfaces. Figures 13 and 14 show 
front views for the total pressure contours on the 
wing surface and in cross flow planes at x = 1.03 
and x = 1.22, and the instantaneous streamlines 
of the inboard tail position. Although, the vortex 
breakdown location is approximately at the same 
position as the uncoupled case, the shape and 
traces of the breakdown flow are different which 
show the upstream effect of the twin-tail motion. 
The vortex cores are moved more upward than 
that of the uncoupled case and continue moving 
upward as the flow travels downstream. Figure 
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15 shows the distribution of the surface pressure 
coefficient covering the wing from x = 0.3 to 
x = 1.0. Figure 16 shows the distribution of the 
leading-edge total structural deflection and the 
root bending moment for the left and right tails 
for 60 dimensionless time after the initial condi- 
tions. The tails deflections and levels of loads 
are higher than those of the uncoupled case. The 
tails deflections are in first, second and third mode 
shapes. 

Figure 17 shows three-dimensional views of 
the leading-edge vortex cores particle traces and 
iso-total pressure surfaces. Figures 18 and 19 show 
front views for the total pressure contours on the 
wing surface and in cross flow planes at x = 1.03 
and x = 1.22, and the instantaneous streamlines 
of the outboard tail position. The shape of the 
vortex breakdown is slightly different than that of 
the uncoupled case. This is because of the lower 
levels of the tail deflections which reduce the up- 
stream effect on the flow. Figure 20 shows the 
distribution of the surface pressure coefficient cov- 
ering the wing from x = 0.3 to x = 1.0. Figure 
21 shows the distribution of the leading-edge to- 
tal structural deflection and the root bending mo- 
ment for the left and right tails for 60 dimension- 
less time after the initial conditions. The tails de- 
flections and levels of loads are higher than those 
of the uncoupled case but still much lower than 
those of the uncoupled and coupled cases of the in- 
board twin-tail position. The tails are deflected in 
one direction only in first and second mode shapes. 

5 CONCLUDING REMARKS 

The effect of coupled and uncoupled bending 
and torsion modes on the twin-tail buffet response 
are investigated for different spanwise positions of 
the twin-tail. The coupled bending and torsion 
modes produce higher deflections and loads than 
those of the uncoupled modes cases. The inboard 
position of the twin-tail produces the largest bend- 
ing and torsion loads and deflections when com- 
pared with the results of the outboard position. It 
has been shown that the larger the tail deflections 
are, the higher the upstream effect on the vortex 
breakdown flow upstream of the tails. The compu- 
tational results presented are in good quantative 
agreement with the experimental data of Wash- 
burn, et.al. (1993). 
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Figure 1: Three-dimensional grid topology and blow-up of the twin tail-delta wing configuration (the 
tails are in rnidspan position). 



Figure 2: Three-dimensional views showing the total pressure on the surfaces, vortex core particle traces 
and iso-total pressure surfaces. Uncoupled case after it = 9,600, Inboard position. 



(a) x = 1.03 (b) x = 1.22 

Figure 3: Snap shots of total pressure in a cross-flow planes. Uncoupled case after it = 9,600, Inboard 
position. 
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Figure 7: Three-dimensional views showing the total pressure on the surfaces, vortex core particle traces 
and iso-total pressure surfaces. Uncoupled case after it = 9,600, Outboard position. 
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Figure 8: Snap shots of total pressure in a cross-flow planes. Uncoupled case after it = 9, 600, Outboard 


Figure 9: Snap shots of instantaneous streamlines in a cross-flow planes. Uncoupled case after it 
Outboard position. 







Figure 10: Distribution of Coefficient of pressure. Uncoupled case after it = 9,600, Outboard position. 



Figure 11: Tail leading-edge total structural deflections and root bending moment for an uncoupled 
bending- torsion case. M ^ = 0.3, a = 30°, R e = 1.25a:10 6 , Outboard position. 



Figure 12: Three-dimensional views showing the total pressure on the surfaces, vortex core particle 
traces and iso-total pressure surfaces. Coupled case after it = 9,600, xg = 0.003, Inboard position. 
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Figure 16: Tail leading-edge total structural deflections and root bending moment for coupled bending- 
torsion case, xg = 0.003. M ^ = 0.3, a = 30°, R e = 1.25a:10 6 , Inboard position. 



Figure 17: Three-dimensional views showing the total pressure on the surfaces, vortex core particle 
traces and iso-total pressure surfaces. Coupled case after it = 9,600, xg — 0.003, Outboard position. 



(a) x = 1.03 (b) x = 1.22 

Figure 18: Snap shots of total pressure in a cross-flow planes. Coupled case after it — 9,600, xg = 0.003, 
Outboard position. 
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Figure 19: Snap shots of instantaneous streamlines in a cross-flow planes. Coupled case after it — 9,600, 
xq = 0.003, Outboard position. 
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Figure 20: Distribution of Coefficient of pressure. Coupled case after it = 9, 600, xq — 0.003, Outboard 
position. 
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ABSTRACT 

The buffet response of the flexible twin- 
tail/delta wing configuration-a multidisciplinary 
problem is solved using three sets of equations on a 
multi-block grid structure. The first set is the un- 
steady, compressible, Reynolds-averaged Navier- 
Stokes equations which are used for obtaining 
the flow-filed vector and the aerodynamic loads 
on the twin tails. The second set is the coupled 
aeroelastic equations which are used for obtain- 
ing the bending and torsional deflections of the 
twin tails. The third set is the grid-displacement 
equations which are used for updating the grid co- 
ordinates due to the tail deflections. The configu- 
ration is pitched at wide range of angles of attack; 
15° to 40°, and the freestream Mach number and 
Reynolds number are 0.3 and 1.25 million, respec- 
tively. With the twin tails fixed as rigid surfaces, 
the problem is solved for the initial flow condi- 
tions. Next, the problem is solved for the twin 
tail response for uncoupled bending and torsional 
vibrations due to the unsteady loads produced 
by the vortex breakdown flow of the leading-edge 
vortex cores. The configuration is investigated for 
two spanwise positions of the twin tails; inboard 
and outboard locations. The computational re- 
sults are validated and are in very good agreement 
with the experimental data of Washburn, et. al. 


INTRODUCTION 

In order to maximize the effectiveness of the 
fighter aircraft that operate well beyond the buffet 
onset boundary, the design of the new generation 
of fighter aircraft should account for both high 
maneuver capabilities at high and wide range of 


angles of attack, and the aeroelastic buffet charac- 
teristics at high alpha. The maneuver capabilities 
are achieved, for example in the F/A-18 fighter, 
through the combination of the leading-edge ex- 
tension (LEX) with a delta wing and the use of 
vertical tails. The LEX maintains lift at high 
angles of attack by generating a pair of vortices 
that trail aft over the top of the aircraft. The vor- 
tex entrains air over the vertical tails to maintain 
stability of the aircraft. At some flight condi- 
tions, the vortices emanating from the highly- 
swept LEX of the delta wing breakdown before 
reaching the vertical tails which get bathed in a 
wake of unsteady highly-turbulent, swirling flow. 
The vortex-breakdown flow produces unsteady, 
unbalanced loads on the vertical tails and causes 
a peak in the pressure spectrum that may be 
tuned to different structural modes depending on 
the angle of attack and dynamic pressure. This 
in turn produce severe buffet on the tails and has 
led to their premature fatigue failure. Therefore, 
the evaluation of the buffet characteristics must 
account for the turbulent characteristics of the 
oncoming flow. If the power spectrum of the tur- 
bulence is accurately predicted, the intensity of 
the buffeting motion can be computed and the 
structural components of the aircraft can be de- 
signed accordingly. 

Experimental investigation of the vertical 
tail buffet of the F/A-18 models have been con- 
ducted by several investigators such as Sellers, et 
al. 1 , Erickson, et al. 2 , Wentz 3 , Lee and Brown 4 , 
and Cole, et al. 5 . These experiments showed that 
the vortex produced by the LEX of the wing 
breaks down ahead of the vertical tails at an- 
gles of attack of 25° and higher producing un- 
steady loads on the vertical tails, and the buffet 
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response occurs in the first bending mode, in- 
creases with increasing dynamic pressure and is 
larger at M = 0.3 than that at a higher Mach 
numbers. Bean and Lee 6 showed that buffeting 
in the torsional mode occurred at a low'er angle 
of attack and at larger levels compared to the 
fundamental bending mode. 

An extensive experimental investigation has 
been conducted to study vortex-tail interaction 
on a 76^ sharp-edged delta wung with vertical 
twin-tail configuration by Washburn, Jenkins and 
Ferman 7 . The vertical tails were placed at nine 
locations behind the wing. The experimental 
data showed that the aerodynamic loads are more 
sensitive to the chordwise tail location than its 
spanwise location. As the tails were moved later- 
ally toward the vortex core, the buffeting response 
and excitation were reduced. Although the tail 
location did not affect the vortex core trajectories, 
it affected the location of vortex-core breakdown. 
Moreover, the investigation showed that the pres- 
ence of a flexible tail can affect the unsteady 
pressures on the rigid tail on the opposite side of 
the model. 

Kandil, Kandil and Massey 8 presented the 
first successful computational simulation of the 
vertical tail buffet using a delta wing-single flex- 
ible vertical tail configuration. The tail was al- 
lowed to oscillate in bending modes. The flow 
conditions and wing angle of attack have been se- 
lected to produce an unsteady vortex- breakdown 
flow. Unsteady vortex breakdown of leading-edge 
vortex cores w r as captured, and unsteady pressure 
forces were obtained on the tail. Kandil, et al. 9 ” 11 
allowed the vertical tail to oscillate in both bend- 
ing and torsional modes. The total deflections 
and the frequencies of deflections and loads of the 
coupled bending-torsion case were found to be 
one order of magnitude higher than those of the 
bending case only. The loads on the tail in the 
transonic flow regime were one order of magni- 
tude lower than those of the subsonic flow. Also, 
it has been shown that the tail oscillations change 
the vortex breakdown locations and the unsteady 
aerodynamic loads on the wing and tail. 

In recent papers by the present authors 12,13 , 
the buffet response of the F/A-18 and a generic 
F-117 tw r in tails w 7 ere considered at a = 30° and 
for different spanwise locations of the twin tails. 
A multi- block grid structure was used to solve 


the problem. The loads, deflections, frequencies 
and root bending moments w^ere reduced as the 
twun tails moved laterally toward the vortex core. 
The outboard location of the tails produced the 
least of these responses. The computational re- 
sults were in full qualitative agreement with the 
experimental data of Washburn, et al. 7 . 

In this paper, we consider the buffet re- 
sponse of delta- wing/twin-tail configuration sim- 
ilar to the one used by Washburn, et. al. 7 . The 
Baldwin and Lomax two-layer turbulent algebraic 
model 14 is used to model flow 7 turbulence. A 
multi-block grid structure is used to solve the 
problem over a wude range of angles of attack 
from 15° to 40°, and for two spanwise locations 
of the twin tails. The computational results are 
compared with the experimental data of Wash- 
burn, et. al. 


FORMULATION 

The formulation consists of three sets of 
governing equations along with certain initial and 
boundary conditions. The first set is the un- 
steady, compressible, Reynolds-averaged Navier- 
Stokes equations. The second set consists of the 
aeroelastic equations for bending and torsional 
modes. The third set consists of equations for 
deforming the grid according to the twin tail de- 
flections. Next, the governing equations of each 
set along with the initial and boundary conditions 
are given. 

Fluid-Flow Equations: 

The conservative form of the dimensionless, 
unsteady, compressible, full Navier-Stokes equa- 
tions in terms of time-dependent, body-conformed 
coordinates £ 2 and £ 3 is given by 

^ 1-3,3 = 1-3 (1) 

dt dt m d£ 3 

where 

* m =* m (Xi,I 2f I 3 ,0 (2) 

Q = j[p,pui,pu 2 ,pu 3 ,pe] t , (3) 

E m and ( E v ) 3 are the £ m -inviscid flux and 
^-viscous and heat conduction flux, respectively. 
Details of these fluxes are given in Ref. 8. The 
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0 ) 


details of the two-layer turbulent algebraic model 
are given in Ref. 14. 

Aeroelastic Equations: 

The dimensionless, linearized governing equa- 
tions for the coupled bending and torsional vibra- 
tions of a vertical tail that is treated as a can- 
tilevered beam are considered. The tail bending 
and torsional deflections occur about an elastic 
axis that is displaced from the inertial axis. These 
equations for the bending deflection, w, and the 
twist angle, 8 , are given by 


d 2 

dz 1 


d 2 e , 


. ,d 2 w 


+m{z)x 9 (z)^(z,t) = N{zJ) (4) 


0 

dz 


c ,«4 

d 2 6 , 


- m(z)x e {z)^(z,t) 




(5) 


where z is the vertical distance from the 
fixed support along the tail length, /*, El and GJ 
the bending and torsional stiffness of the tail sec- 
tion, m the mass per unit length. Is the mass- 
moment of inertia per unit length about the elastic 
axis, x$ the distance between the elastic axis and 
inertia axis, N the normal force per unit length 
and M t the twisting moment per unit length. 
The characteristic parameters for the dimension- 
less equations are c*, a^, p ^ and c* ja*^ for the 
length, speed, density and time; where c* is the 
delta wing root-chord length, a ^ the freestream 
speed of sound and p ^ the freestream air density. 
The geometrical and natural boundary conditions 
on w and 6 are given by 


by 


, . . dw , „ , d 2 w /1 

u?(0, 0 = - 7 ^— (0,0= u 0 


d_ 

dz 

88, 




= 0 ( 6 ) 
(U 


0(0, t) = Q;(lut) = 0 
The solution of Eqs. (4) and (5) are given 


w[z,t) = ^r<t>,(z)q,(t) 


( 8 ) 


i = l 


XT 

*O «0 = <Pj(:)qj(t) 

i=/+i 

where and d> 3 are comparison functions 
satisfying the free- vibration modes of bending and 
torsion, respectively, and q t and qj are generalized 
coordinates for bending and torsion, respectively. 
In this paper, the number of bending modes, /, 
is six and the number of torsion modes, A/ - /, 
is also six. Substituting Eqs. (8) and (9) into 
Eqs. (4) and (5) and using the Galerkin method 
along with integration by parts and the boundary 
conditions, Eqs (6) and (7), we get the following 
equation for the generalized coordinates q t and qj 
in matrix form: 


Mu A / 12 / ix \ , A'u 0 ( q t \ 

A /21 A /22 . \ / . 0 A '22 \ qj ) 



where 

Mu = /o' m<!> r <pidz 

M\2 = M2 1 — Jo mx s 0 r Ojdz > (II) 

4/22 — Jq I$<f>3<Pjdz _ 


K n = & El^fydx \ 
K 22 = j‘‘GJ d -t^dz j 

til = Jo* 0 rNdz 1 
N 2 = J l o' d>,M t dz j 


( 12 ) 

(13) 


Similar aeroelastic equations were devel- 
oped for sonic analysis of wing flutter by 
Strganac 15 . The numerical integration of Eqs. 
(ll)-(13)is obtained using the trapezoidal method 
with 125 points to improve the accuracy of inte- 
grations. The solution of Eq. (10), for q t ;i = 
1,2,....,/, and q 3 \j = /+ 1,...., A/, is obtained us- 
ing the Runge-Kutta scheme. Next, w, and 8 are 
obtained from Eqs. (8) and (9). 


Grid Displacement Equations: 

Once w and 8 are obtained at the n + 1 time 
step, the new grid coordinates are obtained us- 
ing simple interpolation equations. In these equa- 
tions, the twin tail bending displacements, 
and their displacement through the torsion angle. 
8^Jl are interpolated through cosine functions. 
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Boundary and Initial Conditions: 

Boundary conditions consist of conditions for 
the fluid flow and conditions for the aeroelastic 
bending and torsional deflections of the twin tail. 
For the fluid flow, the Riemann-invariant bound- 
ary conditions are enforced at the inflow and out- 
flow boundaries of the computational domain. At 
the plane of geometric symmetry, periodic bound- 
ary conditions are specified with the exception of 
grid points on the tail. On the wing surface, the 
no-slip and no-penetration conditions are enforced 
and = 0. On the tail surface, the no-slip and 
no-penetration conditions for the relative velocity 
components are enforced (points on the tail sur- 
face are moving). The normal pressure gradient is 
no longer equal to zero due to the acceleration of 
the grid points on the tail surface. This equation 
becomes j£ = — pa t .h, where a t is the accelera- 
tion of a point on the tail and ft is the unit normal. 

Initial conditions consist of conditions for 
the fluid flow and conditions for the aeroelastic 
deflections of the twin tail. For the fluid flow, 
the initial conditions correspond to the freestream 
conditions with no-slip and no-penetration condi- 
tions on the wing and tail. For the aeroelastic 
deflections of the tail, the initial conditions for 
any point on the tail are that the displacement 
and velocity are zero, u;(^,0) — 0, ^7(2,0) = 0, 
0{z, 0) = 0 and Jf(*,0) = 0. 


METHOD OF SOLUTION 

The first step is to solve for the fluid flow 
problem using the vortex- breakdown conditions 
and keeping the tail as a rigid beam. Navier- 
Stokes equations are solved using the implicit, 
flux-difference splitting finite- volume scheme. The 
grid speed is set equal to zero in this step. 
This step provides the flow field solution along 
with the pressure differences across the tails. The 
pressure differences are used to generate the nor- 
mal force and twisting moment per unit length 
of each tail. Next, the aeroelastic equations are 
used to obtain the twin tail deflections, Wij % k an< i 
0 t ■ The grid displacement equations are then 
used to compute the new grid coordinates. The 
metric coefficient of the coordinate Jacobian ma- 
trix are updated as well as the grid speed, 

This computational cycle is repeated every time 
step. 


COMPUTATIONAL APPLICATIONS 
AND DISCUSSION 

Twin Tail-Delta Wing Configuration: 

The twin tail-delta wing configuration con- 
sists of a 76°-swept back, sharp-edged delta wing 
(aspect ratio of one) and dynamically scaled flex- 
ible twin tails similar to those used by Washburn, 
et. al. 7 . The vertical tails are oriented normal 
to the upper surface of the delta wing and have 
a centerline sweep of 53.5°. A multi-block grid 
consisting of 4 blocks is used for the solution of 
the problem. The first block is a O-H grid for the 
wing and upstream region, with 101X50X54 grid 
points in the wrap around, normal and axial di- 
rections, respectively. The second block is a H-H 
grid for the inboard region of the twin tails, with 
23X50X13 grid points in the wrap around, nor- 
mal and axial directions, respectively. The third 
block is a H-H grid for the outboard region of the 
twin tails, with 79X50X13 grid points in the wrap 
around, normal and axial directions, respectively. 
The fourth block is a O-H grid for the down- 
stream region of the twin tails, with 101X50X25 
grid points in the wrap around, normal and axial 
directions, respectively. Figure 1 shows the grid 
topology and a front view blow-up of the twin 
tail-delta wing configuration. 

Each tail is made of a single Aluminum spar 
and Balsa wood covering. The Aluminum spar has 
a taper ratio of 0.3 and a constant thickness of 
0.001736. The chord length at the root is 0.03889 
and at the tip is 0.011667, with a span length 
of 0.2223. The Aluminum spar is constructed 
from 606 1-T6 alloy with density, p, modulii of 
elasticity and rigidity, E and G, of 2693 kg/m 3 , 
6.896X10 10 iV/m 2 and 2.5925A"10 l ° N/m 2 \ re- 
spectively. The corresponding dimensionless 
quantities are 2198, 4.595 x 10 5 and 1.727 x 10 5 ; 
respectively. 

The Balsa wood covering has a taper ra- 
tio of 0.23 and aspect ratio of 1.4. The chord 
length at the root is 0.2527 and at the tip is 
0.058, with a span length of 0.2223. The Balsa 
thickness decreases gradually from 0.0211 at the 
tail root to 0.0111 at the tail midspan and then 
constant thickness of 0.0111 is maintained to 
the tail tip. The tail cross section is a semi- 
diamond shape with bevel angle of 20°. The 
Balsa density, modulii of elasticity and rigidity, E 
and G, are \ 79.7 kg/m 3 , 6.896A'10 8 N/m 2 and 
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2.5925AT0 8 Njm 2 ; respectively. The correspond- 
ing dimensionless quantities are 147, 4.595 x 10 3 
and 1.727 x 10 3 ; respectively. Figure 2 shows a 
schematic of the tail used in this study. The tails 
are assumed to be magnetically suspended and 
the leading edge of the tail root is positioned at 
x/c — 1.0, measured from the wing apex. The 
configuration is pitched at a wide range of angles 
of attack 15° to 40°, and the freestream Mach 
number and Reynolds number are 0.3 and 1.25 x 
10 6 ; respectively. 

The configuration is investigated for two 
spanwise positions of the twin tails; the inboard 
location and the outboard location corresponding 
to a separation distance between the twin tails of 
33% and 78% of the wing span; respectively. 


RESULTS AND DISCUSSION 
Parametric Study: 

Figure 3 shows the effect of the spanwise 
tail location on the configuration lift coefficient as 
a function of angle of attack. The lift coefficient 
is almost insensitive to the spanwise tail location. 
The figure also shows a very close agreement with 
the experimental results. Figures 4a and 4b show 
the mean root bending moments and the RMS 
root bending moments as a function of angle of 
attack, respectively. Positive moments correspond 
to an outward force on the tails. The mean value 
of the moment increases with the angle of attack 
and then reduces at high angles. The RMS value 
of the moment also increases with the angle of 
attack. The buffet loads in terms of the root 
bending moment are greatest with the tails in the 
inboard position; almost throughout the whole 
range of angle of attack. The agreement with 
the experimental data is good. The discrepancies 
in the results are attributed to the fact that the 
structural model of the tail is not fully identical to 
the experimental model. The additional Ballast 
weights in the tails of the experimental model is 
not modeled, and Washburn in his experimental 
work used one flexible tail and one rigid tail while 
we are using two flexible tails. In the experimen- 
tal work, the presence of a flexible tail was found 
to affect the loads and pressures on the other rigid 
tail. 

Figure 5 shows the normalized RMS of the 
surface pressures (p/q,^) on the inner and outer 


surfaces of the inboard tails at the specified five 
locations shown in Fig. 2. The experimental re- 
sults of Washburn 7 are also shown. The RMS 
surface pressures in all locations increases with 
the increase of angle of attack, where the vortex 
breakdown moves upstream of the tails. The outer 
surface of the tails experience larger RMS pres- 
sure levels than those of the inner surface. On 
the inner surface, the nearly tip point (location 1) 
experiences larger pressure levels than the nearly 
root point (location 5). On the outer surface, the 
nearly root point experiences larger pressure lev- 
els than those of the nearly tip point. Figure 6 
shows the normalized RMS of the surface pres- 
sures on the inner and outer surfaces of the out- 
board tails at the specified five locations. The ob- 
servations are similar to the case of inboard tails, 
except that the inner surface RMS pressure levels 
are larger than those of the outer surface levels 
of the outboard tail case. The inboard tails ex- 
perience larger pressure levels than those of the 
outboard tails. These results are in good qualita- 
tive agreement with the experimental data. The 
discrepancies between the two results are, again, 
attributed to the reasons stated before. 

Detailed Results at a = 25°: 

Figures 7 and 8 show front views for the 
total pressure contours on the wing surface and 
in cross flow planes at x = 1.03 and x = 1.22, and 
the instantaneous streamlines of the inboard tail 
position. The vortex cores are almost symmetric 
at both locations, and they are totally outboard 
of the twin tail. The cores are moved upward as 
the flow traveled downstream. Smaller size vortex 
cores appear underneath the primary wing vortex 
and it becomes larger in size as it travels down- 
stream. These are the tail vortices observed by 
Washburn 7 . The tail vortices exist at the outer 
surfaces of the tails and they are rotating in the 
opposite direction to those of the primary wing 
vortices. Figure 9 shows the distribution of the 
total structural deflection and the root bending 
moment for the left and right tails. The tails 
deflections are in first, second and third mode 
shapes. Both deflections and root bending mo- 
ment increase rapidly and reach periodic response 
around V* — 20. 

Figures 10 and 11 show front views for the 
total pressure contours on the wing surface and 
in cross flow planes at x = 1.03 and x — 1.22, and 
the instantaneous streamlines of the outboard tail 
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position. The tails cut through the vortex break- 
down flow of the leading-edge vortex cores. The 
flow is almost symmetric. The tail vortices are 
also outboard of the tails but larger in size than 
those of the inboard tails case. The location of 
the vortex core with respect to the tail produces 
an increase in the aerodynamic damping, causing 
the tail deflection to decrease. The tail vortices 
are also shown to rotate in the opposite direction 
to those of the primary wing vortices. Figure 12 
shows the distribution of the total structural de- 
flection and the root bending moment for the left 
and right tails. The levels of deflection are lower 
than those of the inboard tails position. The fre- 
quency of the root bending moment is higher than 
that of the inboard tail case, and the deflections 
seem to be damped and stable at this angle of 
attack. Figure 13 shows the distribution of the 
surface pressure coefficient covering the wing from 
x = 0.3 to x = 1.0. Typical turbulent flow dis- 
tribution are observed, where the largest suction 
peaks are pronounced at the position of the wing 
vortex cores. It is obvious that the vortex break- 
down does not occur over the wing surface at this 
angle of attack. 


Detailed results at a = 40°: 

Figures 14 and 15 show front views for the 
total pressure contours on the wing surface and 
in cross flow planes at x = 1.03 and x = 1.22, and 
the instantaneous streamlines for the inboard tail 
position. The primary leading-edge vortex cores 
experience a breakdown at about x — 0.28, and 
the vortex breakdown flow becomes large in size, 
and partially covers the region inboard the twin 
tails. The vertical position of the vortex cores is 
higher than that at a = 25°, and it moves upward 
as the flow travels downstream. The tail vortices 
are more distinct than those of a = 25°. Figure 
16 shows the distribution of the total structural 
deflection and the root bending moment for the 
left and right tails. The two tails are deflected 
in one direction only in first, second and third 
mode shapes. The frequency of the root bending 
moment is lower than that of the case of a = 25°. 
Periodic responses have not been reached within 
the computational time covered. 

Figures 17 and 18 show front views for the 
total pressure contours on the wing surface and 
in cross flow planes at x = 1.03 and x = 1.22, and 
the instantaneous streamlines for the outboard 


tail position. The tails cut through the vortex 
breakdown flow of the leading-edge vortex cores, 
and the breakdown is larger in size than that of 
the case of a = 25°. The vortex cores are moved 
more upward than that of the case of a — 25° 
and continue moving upward as the flow travels 
downstream. Also, the tail vortex is larger in 
size than that of the a = 25° case. This also 
would increase the aerodynamic damping on the 
tails. The vortex breakdown flow is almost sym- 
metric. Figure 19 shows the distribution of the 
total structural deflection and the root bending 
moment for the left and right tails. The two tails 
are deflected in one direction only in first and 
second mode shapes. The frequency of deflection 
is the same as that of the inboard tail case but 
lower than that of the ca se of a = 25°. The load 
levels are twice as high as those of the case of 
a = 25°. Figure 20 shows the distribution of 
the surface pressure coefficient covering the wing 
from x — 0.3 to x = 1.0. Suction peaks observed 
over the wing are higher than those of the case 
of a = 25°. It is indicated from the distribution 
that vortex breakdown covers most of the wing. 
After 9,600 time steps, the distribution changed 
completely due to the upstream effect of the tail 
deflections, and the vortex breakdown covers the 
whole wing. 

Figures 21 and 22 show the tail surface flow 
on the inner and outer surfaces of the outboard 
tails, respectively, for a = 20°, 30° and 40°. The 
flow separation line is shown to start from the 
leading edge of the root to the middle of the tail 
tip. The separation line moves downstream as the 
angle of attack increases. The tail vortex attach- 
ment line is observed on the outer surface of the 
tail near and parallel to the tail root and along 
the tail root. This is more clear at higher angles 
of attack, a = 30° and 40°. These results are in 
complete agreement with the experimental data 
of Washburn 7 . 

CONCLUDING REMARKS 

The computational results of the FTNS3D 
code of the present paper and the experimental 
data of Washburn, et. al. are in very good agree- 
ment. In general, the RMS levels of the loads 
increase as the angle of attack increases. How- 
ever, the case of inboard tails at a — 25° pro- 
duces the largest buffet loads and deflections. The 
inboard location of the twin tails produces the 
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largest bending-torsion loads, deflections, frequen- 
cies and root bending moments when compared 
with the midspan and outboard locations. The 
frequencies of the loads and deflections of the tails 
decrease as the angle of attack increases. 
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Figure 1: Three-dimensional grid topology and blow-up of the twin tail-delta wing configuration 
(the tails are in midspan position). 




Figure 3: Effect of spanwise tail location 
on the lift coefficient curve. 




Figure 4: Mean and RMS root bending moment on the flexible tails. 
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RMS pressures, inboard tails, P/q RMS pressures, inboard tails, P/q 




(a) Inner surface (b) Outer surface 


Figure 5: RMS pressures at the specified five locations on the tail, inboard position. 




Figure 6: RMS pressures at the specified five locations on the tail, outboard position. 
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Figure 7: Snap shots of total pressure in a cross-flow planes, (Inboard position). = 0.3, a = 25", 
R t = 1 .25x10®. 
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(a) Left tail 


(b) Right tail 


Figure 19: Total structural deflections and root bending moment for an uncoupled bending-torsion case. 
A/co = 0.3, a = 40°, R e = 1.25xl0 6 , (Outboard position). 
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(a) Initial conditions 


(b) Uncoupled case, after it 


Figure 20: Distribution of Coefficient of pressure. Outboard position, Mrc — 0.3, < 
1.25il0 6 . 


(a) a = 20° 


(b) a = 30° 


(c) a = 40° 


Figure 21: Tail surface flow, inner surface. Uncoupled case after it = 9,600, outboard tails 





(a) a = 20° 


(b) a = 30 fl 


(c) a = 40° 


Figure 22: Tail surface flow, outer surface. Uncoupled case after it = 9,600, outboard tails 
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Abstract 

Aeronautical problems and applications of multidisci- 
plinary, time-dependent nature are addressed in this paper. 
Formulations of these problems, which include governing 
equations, boundary and initial conditions, are presented. For 
fiuids/stmctures interaction problems, three sets of governing 
equations are used. These are the unsteady, Navier-Stokes 
(NS) equations, the aeroelastic equations of coupled or uncou- 
pled bending and torsion equations and the grid displacement 
equations. The latter equations are used to move the compu- 
tational grid due to the structural elastic deformations. For 
fluids/dynamics interaction problems, two sets of governing 
equations are used. These are the unsteady NS equations and 
the Euler equations of rigid-body dynamics. These sets of 
equations are sequentially solved time accurately using the up- 
wind, flux-difference splitting scheme for the fluids equations 
and the four-stage Rungc-Kutta scheme for the aeroelastic 
and dynamics equations. The grid displacement equations are 
algebraic equations which are solved to obtain the new grid 
coordinates. 

Computational applications and validations cover delta 
wings/ twin tail configurations simulating the tail buffet prob- 
lem and delta wings undergoing natural mode response in 
rolling oscillations simulating the wing rock problem. Vali- 
dation of the computational results are also presented and dis- 
cussed. 

Introduction 

The first multidisciplinary aeronautical problem is that of 
fluids/structures interaction. Vertical tails buffet represents 
a challenging aeronautical application in this area. Modem 
aircraft are designed to fly and maneuver high angles of attack 
and at high loading conditions. This is achieved, for example 
in the F/A-18 fighter, through the combination of a leading- 
edge extension (LEX) with a delta wing and the use of highly 
swept-back twin vertical tails. The LEX maintains lift at high 
angles of attack by generating a pair of vortices that trail aft 
over the top of the delta wings. The vortex entrains air over 
the vertical tails to maintain stability of the aircraft. This 
combination of LEX, delta wing and vertical tails leads to the 
aircraft excellent agility. However, at some flight conditions, 
the vortices emanating from the highly-swept LEX of the delta 
wing breakdown before reaching the vertical tails, which get 
bathed in a wake of unsteady highly-turbulent, swirling flow. 
The vortex-breakdown flow produces unsteady, unbalanced 
loads on the vertical tails which in turn produce severe buffet 
of the tails. The buffet loads result into premature fatigue 
failure of the tails. 
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Experimental investigation of the vertical tail buffet of the 
F/A-18 models have been conducted by several investigators 
such as Sellers et al (l) ., Erickson et al (2) ., Wentz (3) and Lee and 
Brown (4) . These experiments showed that the vortex produced 
by the LEX of the wing breaks down ahead of the vertical tails 
at angles of attack of 25° and higher and that the breakdown 
flow produced unsteady loads on the vertical tails. Cole, Moss 
and Doggett™ tested a rigid. 1/6 size, full-span model of an 
F-18 airplane that was fitted with flexible vertical tails of two 
different stiffness. Vertical-tail buffet response results were 
obtained over the range of angle of attack from -10* to +40° , 
and over the range of Mach numbers from 0.3 to 0.95. Their 
results indicated that the buffet response occurred in the first 
bending mode, increased with increasing dynamic pressure and 
was larger at M = 0.3 than that at a higher Mach number. 

An extensive experimental investigation has been con- 
ducted to study vortex-twin tail interaction on a 76* sharp- 
edged delta wing with vertical twin-tail configuration by Wash- 
bum, Jenkins and Ferman <6 \ The vertical twin tails were 
placed at nine locations behind the wing. The experimental 
data showed that the aerodynamic loads are more sensitive to 
the chordwise tail location than its spanwise location. As the 
tails were moved toward the vortex core, the buffet response 
and excitation were reduced. Although the tail location did 
not affect the vortex core trajectories, it affected the location 
of vortex -core breakdown. Moreover, the investigation showed 
that the presence of a flexible tail can affect the unsteady pres- 
sures on the rigid tail on the opposite side of the model. In 
a recent study by Bean and Lee™ tests were performed on a 
rigid 6% scale F/A-18 in a trisonic blowdown wind tunnel over 
a range of angle of attack and Mach number. The flight data 
was reduced to a non-dimensional buffet excitation parameter, 
for each primary mode. It was found that buffeting in the tor- 
sional mode occurred at a lower angle of attack and at larger 
levels compared to the fundamental bending mode. 

Tail buffet studies were also conducted on a full-scale, 
production model F/A-18 fighter aircraft in the 80-by— 120 foot 
wind tunnel at NASA Ames Research Center by Meyn and 
James™ and Pettit, Brown and Pendleton™. The test matrix 
covered an angle of attack range of 18° to 50° and a side-slip 
range of -16° to 16° with wind speed up to 100 Knots. The 
maximum speed corresponds to a Reynolds number of 1.23 x 
10 7 and a Mach number of 0.15. 

Kandil, Kandil and Massey (I0) presented the first successful 
computational simulation of the vertical flexible tail buffet us- 
ing a delta wing-vertical tail configuration. A 76° sharp-edged 
delta wing has been used along with a single rectangular ver- 
tical tail which was placed aft the wing along the plane of 
geometric symmetry. The flexible tail was allowed to oscil- 
late in bending modes. The flow conditions and wing angle 
of attack have been selected to produce an unsteady vortex- 
breakdown flow. Unsteady vortex breakdown of leading-edge 
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vortex cores was captured, and unsteady pressure forces were 
obtained on the tail. 

Kandil, Massey and Kandi! (,n extended the technique used 
in Ref. 10 to allow the vertical tail to oscillate in both bending 
and torsional modes. The total deflections and the frequencies 
of deflections and loads of the coupled bending- ton ion case 
were found to be one order of magnitude higher than those 
of the bending case only. Also, it has been shown that the 
tail oscillations change the vortex breakdown locations and the 
unsteady aerodynamic loads on the wing and tail. 

Kandil. Massey and Sheta (l2) studied the effect of coupling 
and uncoupling the bending and torsional modes for a long 
computational time, and the flow Reynolds number on the 
buffet response, of a single rectangular flexible tail. It has been 
shown that the coupled response produced higher deflection 
than that of the uncoupled response. Moreover, the response 
of the coupled case reached periodicity faster than that of the 
uncoupled case. It has also been shown that he deflections of 
the iow-Reynolds number case were substantially lower than 
that of the high Reynolds number case. 

In a recent paper by Kandil, Sheta and Massey< 13 \ the 
buffet response of a single swept-back vertical flexible tail in 
transonic flow at two angles of attack (20°. 28°) has been 
studied. It has been shown that the aerodynamic loads and 
bending- torsion deflections of the tail never reached periodic 
response and that the loads were one order of magnitude lower 
than those of Ref. 12 of the subsonic flow. 

In a very recent paper by Kandil, Sheta and Massey {!4 \ 
the buffet response of the F/A-18 twin tails were considered. 
The configuration consisted of a 76° -swept back, sharp-edged 
delta wing and a trai ling-edge -extension on which the F/A-18 
twin tails were attached as cantilevers. A multi-block grid was 
used to solve the problem for two lateral locations of the twin 
tails; a midspan location and an inboard location. 

The second multidisciplinary aeronautical problem is that 
of fluids/dynamics interaction. Wing rock of highly swept-back 
wings represents another challenging aeronautical application, 
in this area 

One frequently encountered lateral instability which limits 
combat effectiveness for all fighter aircraft is the limit-cycle 
rolling oscillation phenomenon known as wing rock. In mod- 
erate to high angle-of-attack dynamic motion, wing rock is 
driven by strong, concentrated vortices originating from the 
leading edges of highly swept wings. Wing rock can occur at 
subsonic airspeeds at angles of attack in the vicinity of stall, 
and at moderate angles of attack in the transonic regime as a 
result of shock-wave/boundary-layer interactions on the wing. 
Generally, the onset of wing rock can be caused by a number 
of different aerodynamic phenomena and is attributed to a loss 
of stability in the lateral/directional mode. 

To understand the wing rock phenomenon, experimental 
investigations have been carried out on simplified delta-wing 
geometries with a single degree of freedom in roll. By avoid- 
ing the complexity of complete aircraft geometries, research 
can focus on the relevant flow physics. Experimental data typ- 
ically consists of flow visualization, time -dependent forces and 
moments and more recently, time-dependent surface pressure 
data. The time -dependent pressure data provides additional in- 
formation that allows for more detailed understanding of the 


mechanisms involved with wing rock which have yet to be 
fully understood. However, these experimental results are lim- 
ited by the difficulties encountered in taking measurements in 
a dynamic environment. 

Computational fluid dynamics (CFD) plays an important 
role in the design process by providing detailed flowfield in- 
formation at a relatively low cost that is unavailable with ex- 
periment alone. A few computational studies have been ini- 
tiated to simulate the wing rock problem. However, due to 
large amounts of computational lime, most of these studies 
have employed various limiting approximations to reduce the 
computational cost. Inherently, these simplifying flow assump- 
tions restrict the applicability of the solution to steady or invis- 
cid flows. For vortical flows where viscous effects dominate, 
computation based on the unsteady Navier-Stokes equations is 
vital. The Navier-Stokes equations can more accurately model 
flow separations, shock development and motion, and shock- 
boundary-layer interaction as well as vortex breakdown and 
vorticity evolution, convection and shedding. 

In 1981, the phenomena of slender wing rock was first 
observed in experiments performed by Nguyen. Using 

an 80° swept delta wing investigation showed that wing rock 
occurred simultaneously with the appearance of asymmetric 
leading-edge vortices. By 1984, Ericsson (,fi) had shown that 
vortex asymmetry could generate wing rock but growth of 
the amplitude was limited by vortex breakdown. Arena {17) 
conducted a thorough experimental investigation of the natural 
response of a slender wing rock in subsonic flow. He identified 
the envelope of damped and self sustaining motion for an 
80° swept wing and qualitatively compared these results with 
computational results. Continuing investigation of wing rock. 
Ng, ef a/., ul) used a water tunnel to compare forced rolling and 
free-to-roll oscillations of delta wings of various sweep angles 
with static conditions. 

Various experimental attempts to control wing rock have 
also been investigated experimentally. Malcolm, et a/. (m 
demonstrated a wing’s rolling moment can be affected by me- 
chanical or pneumatic manipulation of the strength or location 
of the leading-edge vortices. In 1993, Walton and Katz (2n) ex- 
ploited this idea and applied leading edge control flaps to a 
free-to-roll double-delta wing. In 1994, Ng, et al., m) demon- 
strated passive control of an 80° swept delta wing undergoing 
wing rock by using flow dividers. At angles of attack higher 
than 30", suppression of wing rock was achieved. However, at 
lower angles of attack, the divider actually promoted the phe- 
nomenon. Using asymmetric tangential leading-edge blowing, 
Wong, et a/., (22) demonstrated positive post-stall roll control 
for a delta wing at an angle of attack of 55°. With an ac- 
tive roll feedback control algorithm employing a proponional- 
derivative compensator, wing rock was stopped in less than 
one cycle of the limit-cycle motion. 

As in experimental investigations of forced rolling oscilla- 
tions. the focus of computational studies is to be able to predict 
and ultimately control the phenomenon of wing rock. In 1985, 
using an unsteady vortex-lattice method, Konstadinopoulos, et 
a/. <23) numerically simulated the subsonic experimental work 
performed by Nguyen, et a/. <l5) . They determined that the 
leading-edge vortex system became unstable as the angle of 
attack was increased which caused a loss of roll damping at 
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small angle of roll. Improving the methods for numerical sim- 
ulation, in 1989, Nayfeh, et o// 24) proceeded to construct phase 
planes which revealed the general global nature of wing rock by 
discussing stable limit cycles, unstable foci, and saddle points. 
This demonstrated the locations of equilibrium positions. By 
1994, Chaderijian and Schiff etaf 2 ^ had solved for flow over 
a 65° swept delta wing at 30° angle of attack and Mach of 0.27 
that was both forced and free to roll under the influence of the 
instantaneous aerodynamic rolling moment. 


initial and boundary conditions. The first set is the unsteady, 
compressible, full Navier-Stokes equations. The second set 
consists of the aeroelastic equations for bending and torsional 
modes. The third set consists of equations for deforming the 
grid due to the structure deflections. Next, the governing 
equations of each set along with the initial and boundary 
conditions are presented. 

Fluid-Flow Equations 


Numerical simulation for the control of wing rock has been 
performed by various authors primarily, using Euler equations 
assuming locally conical flow. In 1991, after developing the 
Navier-Displacemem equations for grid deformation, Kandil 
and Salman 061 effectively controlled the wing rock response of 
an 80° swept delta wing at 30° angle of attack and Mach num- 
ber of 1.2 by using tuned antisymmetric leading-edge flap oscil- 
lations. They also applied the locally conical Euler equations to 
the same problem at Mach 1 .4. The three-dimensional flow so- 
lution of Euler equations at Mach 0.3 were also considered/ 271 
Noting the loss of aerodynamic damping rolling moment at 
the zero angular velocity value, they determined that the hys- 
terias responses of position and strength of the asymmetric 
right and left primary vortices were responsible for wing rock 
and that the phenomenon could be actively controlled through 
the use of leading edge flaps. In 1991, Kandil and Salman 0 * 1 
solved the thin-layer locally conical Navier-Stokes equations 
for delta wing at 35° angle of attack. It was again showed 
that the wing-rock phenomenon could be controlled by using 
tuned anti-symmetric leading-edge flap oscillations. Details of 
this work were published in Salman’s dissertation/ 291 In 1993, 
Lee-Rausch and Batina 001 also investigated control of wing 
rock using locally conical Euler equations using leading-edge 
flaps. Their study focused on a 75 g swept sharp-edged delta 
wing at a free-stream Mach number of 1.2 at various angles 
of attack. 

In Ref. (31), Menzies and Kandil presented three cases of 
computationally simulated natural rolling response of a delta 
wing in transonic and subsonic flow. This was the only known 
published study in the transonic flow regime using the NS equa- 
tions. Transonic flow over a 65® swept, cropped delta wing 
with breakdown of the leading edge vortices demonstrated self 
sustained rolling oscillations until breakdown dominates the 
flow field. Two cases of subsonic flow over an 80° swept 
wing demonstrated either damped or self-sustained rolling os- 
cillations as a function of angle of attack. A complete inves- 
tigation of the aerodynamic response of the wing, the effects 
of Mach number, angle of attack, and vortex breakdown were 
presented. 

In this paper, the formulation and method of solution of 
each multidisciplinary problem is presented. Computational 
results and validations are also presented and discussed for 
each problem. The computational applications arc focused on 
the high-angle-of-attack aerodynamics of delta wings and delta 
wings/vertical tails configurations. 

Formulations 


Fluid/Structures Interaction 


The conservative form of the dimensionless, unsteady, 
compressible, full Navier-Stokes equations in terms of time- 
dependent, body-conformed coordinates f 2 and £ 3 is given 
by 


8Q dE m 

di 


*(&), 

dz- 


= 0; m = 1 


— 3,5 = 
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( 1 ) 


where 


£"=r(x i,*a,*3.0 


( 2 ) 

(3) 


E m and (E v ) t are the £ m -inviscid flux and £* -viscous and 
heat conduction flux, respectively. Details of these fluxes are 
given in Ref. 10. 


Aeroelastic Equations 

The dimensionless, linearized governing equations for cou- 
pled bending and torsional vibrations of a vertical tail, that is 
treated as a cantilevered beam, are considered. The tail bend- 
ing and torsional deflections occur about an elastic axis that 
is displaced from the inertial axis. These equations for the 
bending deflection, w, and the twist angle, 9, are given by 


d 2 

dz 2 




(4) 


d 7 9 

-/,(*) — ( z,0=-A/e(*,0 

where z is the vertical distance from the fixed support along the 
tail length, /„ El and GJ the bending and torsional stiffness of 
the tail section, m the mass per unit length. Iq the mass-moment 
of inertia per unit length about the elastic axis, x 9 the distance 
between the elastic axis and inertia axis, N the normal force 
per unit length and M t the twisting moment per unit length. 
The characteristic parameters for the dimensionless equations 
are c*, a^, and c*/a^ for the length, speed, density and 
time; where c* is the delta wing root-chord length, a w m the 
freestream speed of sound and p ^ the freestream air density. 
The geometrical and natural boundary conditions on w and 9 
are given by 


-CM) = 




( 6 ) 
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= 0 


For fluids/structures interaction problems, the formulation 
consists of three sets of governing equations along with certain 


*(<U) = ^(U) = ° 


(7) 
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The solution of Eqs. (4) and (5) are given by 

/ 

“(M) = £ 4<(*)«i(0 (8) 

1=1 


*(M) 


M 


i=/+i 


(9) 


where and 4> j are comparison functions satisfying the frce- 
vibration modes of bending and torsion, respectively, and q { 
and q { are generalized coordinates for bending and torsion, 
respectively. In this paper, the number of bending modes, /, 
is six and the number of torsion modes, M - /, is also six. 
Substituting Eqs. (8) and (9) into Eqs. (4) and (5) and using 
the Galerkin method along with integration by parts and the 
boundary conditions, Eqs (6) and (7), we get the following 
equation for the generalized coordinates and q , in matrix 
form: 



(Nx\ ;i = 1*2, / 

\N 2 ) ;j = / + 1,...., A/ 


( 10 ) 


where 


AT, i 


= /o' rn<t> r 4> x dz 

= W.I = /o' mx^^jdz 

= /o' It<P,4>jdz 

*■" = /.'' 1 


i.Md, 


(ID 

( 12 ) 

(13) 


The numerical integration of Eqs. (ll)-(!3) is obtained 
using the trapezoidal method with 125 points to improve the 
accuracy of integrations. The solution of Eq. (10), for t = 

1,2....,/, and j — J+l, A/, is obtained using the Runge- 

Kutta scheme. Next, w, and 9 are obtained from Eqs. (8) and 
(9). 


Grid Displacement Equations 

Once w and 9 are obtained at the n + 1 time step, the new 
grid coordinates are obtained using simple interpolation equa- 
tions. In these equations, the twin tail bending displacements, 
and their displacements through the torsion angle, 9*+l 
are interpolated through cosine functions, with maximum dis- 
placements at the tail surface and zero displacements at the 
computational boundaries. 


(points on the tail surface are moving). The normal pressure 
gradient is no longer equal to zero due to the acceleration of 
the grid points on the tail surface. This equation becomes 

= -pa, ■ n, where <5 t is the acceleration of a point on the 
tail and h is the unit normal. 

Initial conditions consist of conditions for the fluid flow 
and conditions for the aeroelastic deflections of the tails. 
For the fluid flow, the initial conditions correspond to the 
freestream conditions with no-slip and no-penetration condi- 
tions on the wing and tail. For the aeroelastic deflections of 
the tail, the initial conditions for any point on its surface are that 
the displacement and velocity are zero, “'(z.O) = O.TT^.O) - 
O,0(z,O) = 0 and £(z,0) = 0. 

Fluids/Dynamics Interaction 

For fluids/dynamics interaction problems, the formulation 
consists of two sets of governing equations, and initial and 
boundary conditions. The first set is the unsteady, compress- 
ible, full Navier-Stokes equations. This set has already been 
given in Eqs. (l)-(3). The second set is the Euler equations 
of rigid body dynamics. In this paper, we consider only one- 
degree of freedom for the wing motion. For the natural rolling 
response, the resultant external aerodynamic rolling moment 
is equated to the time rate of change of the angular rolling 
moment about the axis of rolling motion. This equation is 
given by 


On,. <( = w, (14) 

where /„ is the principal mass moment of inertia for the wing, 
and w r is the roll accel era lion. 


Boundary and Initial Conditions and Grid Motion 

Since the wing is undergoing rolling motion, the grid is 
moved with the same angular motion as that of the body. 
The grid speed, and the metric coefficient, are 

computed at each time step of the computational scheme. 
Consequently, the kinematic boundary conditions at the inflow- 
outflow boundaries and at the wing surface are expressed in 
terms of the relative velocities. Again, the dynamic boundary 
condition, on the wing surface is no longer equal to zero. 
This condition is modified for the oscillating wing as: 


dp 

8n 


= —pa h 

wing 


05) 


where a is the acceleration of a point on the wing flat surface; 
n, the unit normal to the wing surface. The acceleration is 
given by: 


Boundary and Initial Conditions 

Boundary conditions consist of conditions for the fluid flow 
and conditions for the aeroelastic bending and torsional deflec- 
tions of the tails. For the fluid flow, the Riemann -invariant 
boundary conditions are enforced at the inflow and outflow 
boundaries of the computational domain. At the plane of 
geometric symmetry, periodic boundary conditions is specified 
with the exception of grid points on the tail. On the wing sur- 
face. the no-slip and no- penetration conditions are enforced and 
= 0. On the tail surfaces, the no-slip and no- penetration 
conditions for the relative velocity components are enforced 


a ~ ux r + wr(w x f) 


(16) 


where u? is the angular velocity. Notice that for a rigid body, 
the ‘position vector f, is not a function of time and hence, 
r = r = 0. Finally, the boundary condition for the temperature 
is obtained from the adiabatic boundary condition and is given 
by: 


3T 

dn f u>injj 


(17) 
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Method of Solution 

Fluids/Structures Interaction 

The first step is to solve for the fluid flow problem using 
vortex -breakdown conditions and keeping the tails rigid. The 
NS equations are solved with the grid speed set equal to 
zero. The implicit, upwind, flux -difference splitting, finite- 
volume scheme is used to solve the unsteady, compressible, 
full Navier-Stokes equations. This scheme uses the flux- 
difference splitting of Roe, and a smooth flux limiter is used 
to eliminate oscillation at locations of large flow gradients. 
The viscous and heat flux terms are linearized in time and the 
cross derivative terms are eliminated in the implicit operator 
and retained in the explicit terms. The viscous terms are 
differenced using second-order accurate central differencing. 
The resulting difference equation is approximately factored to 
solve the equations in three sweeps in the and 

directions. The computational scheme is coded in the computer 
program “FTNS3D" This step provides the flow field solution 
along with the pressure difference across each of the tails. 
The pressure difference is used to generate the normal force 
and twisting moment per unit length of each tail. Next, the 
aeroelastic equations, Eqs. (4)-(13), are used to obtain the tails 
deflections, w ijjk and 0,,^*. The grid displacement equations 
are then used to compute the new grid coordinates. Once 
the tails are allowed to deform, the metric coefficients of the 
coordinate Jacobian matrix are updated as well as the grid 
speed, ^2.. This computational cycle is repeated every time 
step. 

Fluids/Dynamics Interaction 

The method of solution consists of three steps. In the first 
step, the problem is solved for the stationary wing at 0° roll 
angle. This solution represents the initial conditions for the 
second step. In the second step, the dynamic initial conditions 
are specified. In this paper, we consider two types of flow 
regime; a transonic case and subsonic cases. For the tran- 
sonic case, the wing is subjected to an initial velocity. For 
the subsonic cases, a quarter cycle of a sinusoidal function is 
specified to roll the wing to a 10* roll angle with zero angular 
velocity while the Navier-Stokes equations are solved accu- 
rately in time. Having specified the dynamic initial conditions, 
the third step proceeds. Applying a four-stage Range- Kutta 
scheme and the specified dynamic initial conditions for 0 and 
0, Eq. (14) is explicitly integrated in time in sequence with 
the fluid dynamic equations. 

Equation (14) is used to solve for 0, 0, and 0 while the 
fluid dynamics equations provide the pressure distribution over 
the wing surface. The pressure distribution is integrated over 
the surface of the wing to determine with respect to 

the axis of geometric symmetry. At each time step, the wing 
and the grid are rotated corresponding to 0 and 0 resulting 
in the natural rolling response of the delta wing to the fluid 
flow. Due to the dynamic nature of the problem, the metric 
coefficients and the grid speed are computed at each time step. 
The computations proceed until periodic response is reached. 

Computational Application And Discussion 
Fluid-Structure Interaction; Twin-Tail Buffet 


The twin tail-delta wing configuration used in this appli- 
cation consists of a 76°-swept back, sharp-edged delta wing 
(aspect ratio of one) and dynamically scaled flexible twin tails 
similar to those used by Washburn, et. al. (Ref 6). The vertical 
tails are oriented normal to the upper surface of the delta wing 
and have a centerline sweep of 53.5° Each tail is made of 
a single Aluminum spar and Balsa wood covering. The Alu- 
minum spar has a taper ratio of 0.3 and a constant thickness of 
0.001736. the chord length at the root is 0.03889 and at the 
tip is 0.01 1667, with a span length of 0.22223. The Aluminum 
spar is constructed from 6061-T6 alloy with density, p, modulii 
of elasticity and rigidity, E and G of 2693 kg/m 3 , 6.896X1 0 ,n 
S/m 2 and 2.5925X10*° N/m 2 \ respectively. The Balsa wood 
covering has a taper ratio of 0.23 and aspect ratio of 1-4. The 
chord length at the root is 0.2527 and at the tip is 0.058, with a 
span length of 0.2223 The Balsa thickness decreases gradually 
from 0.021 1 at the tail root to 0.0111 at the tail midspan and 
then constant thickenss of 0.0 11 1 is maintained to the tail tip. 
The tail cross section is a semi-diamond shape with bevel angle 
of 20°. The Balsa density, modulii of elasticity and rigidity, 
E and G, are 179.7 kg/m 3 , 6.896X10* N/m 2 and 2.5925X1 0 8 
N/m 2 \ respectively. The tails are assumed to be magnetically 
suspended and the leading edge of the tail root is positioned 
at x/c = 1.0, measured from the wing apex. The configura- 
tion is pitched at an angle of attack of 30° and the freestream 
Mach number and Reynolds number are 0.3 and 1.25 x 10 A ; 
respectively. 

A multi-block grid consisting of 4 blocks is used for the 
solution of the problem. The first block is a O-H grid for 
the wing and upstream region, with 101X50X54 grid points 
in the wrap around, normal and axial directions, respectively. 
The second block is a H-H grid for the inboard region of the 
twin tails, with 23X50X13 grid points in the wrap around, 
normal and axial directions, respectively. The third block is 
a H-H grid for the outboard region of the twin tails, with 
79X50X13 grid points in the wrap around, normal and axial 
directions, respectively. The fourth block is a O-H grid for 
the downstream region of the twin tails, with 101X50X25 
grid points in the wrap around, normal and axial directions, 
respectively. Figure I shows the grid topology and a blow-up 
of the twin tail-delta wing configuration. 

The configuration is investigated for three spanwise posi- 
tions of the twin tails; the inboard location, the midspan lo- 
cation and the outbaord location corresponding to a separation 
distance between the twin tails of 33%, 56% and 78% of the 
wing span; respectively. 

Inboard Location of Twin Tails (33% wing span) 

The spanwise distance between the two tails is 33% of the 
wing span. Figure 2 shows three-dimensional and front views 
for the initial conditions with the surface total pressure contours 
and the streamlines of the vortex cores. Figure 3 shows the 
static pressure contours and the instantaneous streamlines in 
a cross flow plane at x = 1.096. The initial conditions are 
obtained after 10,000 time steps, A t - 0.001, with the twin 
tails kept rigid. It is observed that the vortex cores experience 
an almost symmetric breakdown on the wing at about the 
75% chordstation. Downstream of the wing, they are totally 
outside of the space between the twin tails. Smaller size vortex 
cores appear under the vortex breakdown flows and at the 
lower edges of the twin tails. These results exactly match 
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Washburn observations. Figures 4 and 5 show the results for 
the twin tails undergoing uncoupled bending- torsion responses 
after 9,600 time steps from the initial conditions. It is observed 
that the breakdown shapes and locations are affected by the 
twin tail oscillations. The vortex breakdown is now strongly 
asymmetric, and the vortex breakdown flows are still outside 
of the space between the twin tails. These results conclusively 
show the upstream as well as the spanwise effect of the twin 
tail oscillations on the vortex breakdown flow. 

Figures 6-8 show the distribution of deflection and load 
responses along the left and right tails every 2 time units, the 
history of deflection and load responses versus time and the 
total structural deflections and root bending moment for the 
left and right tails. It is observed that the bending and torsion 
responses are in their first and second mode shapes. The 
frequencies of the bending deflections are less than one-half 
those of the torsion deflections. The normal forces are out of 
phase of the bending deflections while the torsion moments are 
in phase with the torsion deflections. The total tail responses 
are in first, second and third mode shapes. Periodic responses 
have not been reached within the computational time covered 
(20,000 time steps =* 10 dimensionless time units). 

Midspan Location of Twin Tails (56% wing span) 

The results of this case are presented in Figs. 9-15. Figures 
9- 1 2 show that the tails cut through the vortex breakdown of the 
leading-edge vortex cores, which are also asymmetric. Figure 
13 shows that the bending deflections are lower than those of 
the inboard case while the torsional deflections are substantially 
lower than those of the inboard case. Moreover, Fig. 13 shows 
that the bending and torsion deflections have a single sign for 
the left and right tails (all are positive or all arc negative). 
Figure 14 shows that both bending and torsion deflections are 
out of phase of the normal force and twisting moment loads. 
The total deflections of Fg. 15 show the same trend. The 
root bending moments of Fig. 15 are also lower than those of 
the inboard case. 

Outboard Location of Twin Ttils (78% wing span) 

Figures 16-20 show the results of this case. Figures 16 
and 17 show that the space between the twin tails include 
larger portion of the vortex breakdown flow of the leading- 
edge vortex cores, than that of the midspan case. The vortex 
breakdown flow is also asymmetric. The vortical flow on the 
lower outside surfaces of the twin tails is larger than any of 
the above two cases. Figures 18-20 show that the bending 
and torsion deflections are lower than those of the midspan 
case. They also show that both bending and torsion deflections 
are out of phase of the bending and torsional loads. The 
frequencies of the bending deflections are still smaller than 
those of the torsion deflections. All these observations are in 
very good agreement with those of Washburn, et. al. (Ref 
6). Figures 21-23 show the histories of the lift and drag 
coefficients versus time for the inboard, midspan and outboard 
locations, respectively. It is observed that the loss in CL is the 
largest for the inboard location case. 

Table !. shows a comparison of the present results of the 
mean root bending moment for flexible twin tails and the lift 


coefficient with rigid twin tails with those of Washburn, et. al. 
(Ref 6), experimental data. 


[ Parameter 

Position 

FTNS3D 

(present) 

WASHBURN (6) 

I Mean Root Bending 
I Moment 

With Flexible Tails 

Inboard 

5.62 X 10 J 

7 43 X 10' 5 

Midspan 

4.22 X 10’ 3 

6.05 X Ifr 5 

Outboard 

3.62 X I0 5 

5 70 X 10 5 | 

flu ft Coefficient 
With Rigid Tails 

Inboard 

1.0423 

1.17 

Midspan 

1.0515 

1.12 

Outboard 

1.0674 

1.17 


Fluids/Dynamics Interaction 

Case I-Transonic Flow over a Cropped Delta Wing 

A 65® swept-back. sharp edged, cropped delta wing of zero 
thickness is considered for the transonic flow solutions. The 
cropping ratio (tip length/root-chord length) is 0.15. An O-H 
grid of 65 x 43 x 105 in the wrap-around, normal, and axial 
directions, respectively, is used. The computational domain 
extends two chord lengths forward and five chord lengths 
downstream of the wing trailing edge. The radius of the 
computational domain is four chord lengths. The minimum 
grid size in the normal direction to the wing surface is 5 x 10^ 
from the leading edge to the plane of symmetry. The initial 
conditions correspond to the solution of the wing held at 20° 
angle of attack and 0° roll angle after 18,000 time steps at a 
Mach number and Reynolds number of 0.85 and 3.23 x 1 0 6 , 
respectively. 

Plots of the initial condition depict a solution characterized 
by weak oblique shocks beneath the primary vortices and a 
strong, transverse terminating shock located a approximately 
x » 0.86 (See Fig. 24). These shocks bound a substantial 
supersonic pocket. Outboard of the oblique shocks, a subsonic, 
separated region depicts a secondary vortex which exists until 
x = 0.91. The primary vortex interacts with the terminating 
shock and enlarges indicating vortex breakdown. The plots 
of the Mach number contours, instantaneous streamlines, and 
surface of constant entropy shown in Fig. 24 depict clearly 
a bubble type vortex breakdown and the flow appears to be 
completely symmetric. 

From the initial conditions, the wing is given an instanta- 
neous roll velocity of 9 = +0.925 x I0 - *. For convention, a 
positive roll velocity indicate that the right hand side of the 
wing when looking in the upstream direction is moving up- 
ward, With a dimensionless mass moment of inertia of l u = 
2.88 x 10" 3 , the wing is free to respond to the rolling mo- 
ment induced by the fluid flow. Figure 25 shows the time his- 
tory of the resultant motion and lift coefficient curve. While 
the motion appears somewhat periodic after t = 30, there is 
a chaotic behavior in the cyclic response due to the vortex 
breakdown which leads to divergence of the motion after five 
cycles of rolling. The highly unsteady nature of the shock in- 
duced vortex breakdown promotes very irregular motion and 
aerodynamic response histories. The lift coefficient indicates 
an initial loss after the onset of motion and fluctuates between 
0.36 and 0.40 during the quasi -periodic response. After t = 
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120 when the wing motion diverges to approximately 24°, the 
lift drops considerably. 

Figure 26 shows snapshots of the Mach contours near the 
wing surface and surfaces of constant entropy depicting the 
primary vortex core and breakdown. During the quasi-periodic 
response, the terminating shock and vortex breakdown location 
oscillate laterally as shown in views a) to d). In contrast 
to subsonic flow where the breakdown and restructuring of 
the vortex serves to dampen the motion;* 171 in transonic flow, 
the sustained oscillation of the wing appears to be caused by 
the asymmetric motion of the breakdown location. However, 
as time progresses, the frequency of oscillation is slightly 
increased. In Ref. 32, it was shown that when forced to 
roll at a reduced frequency of 2x\ the transverse shock was 
weakened as a result of the motion and the breakdown washes 
downstream. In this case, after t = 120, (shown in view e), the 
shock on the upward moving side appears to weaken as a result 
of the increased rolling frequency and the breakdown washes 
downstream. On the downward moving side, the breakdown 
location advances to apex of the wing which results in a 
significant drop in the lift. Without breakdown on the right 
side, the wing rolls until equilibrium is reached at a positive 
mean roll angle of 24.2°. Small amplitude fluctuations of 
this roll angle are observed due to the unsteady nature of the 
complete breakdown on the left side of the wing. 

Case II-Subsonic Flow over a Delta Wing 
undergoing Damped Rolling Oscillations 

In order to compare with available experimental data, an 
80° swept-back, sharp-edged delta wing of zero thickness is 
considered for the subsonic flow solutions. This wing was 
modeled after the experimental model used by Arena.* 171 An 
O-H grid of 65 x 43 x 84 in the wrap-around, normal, and axial 
directions, respectively, is used. The computational domain 
extends two chord lengths forward and five chord lengths 
backward from the wing trailing edge. The radius of the 
computational domain is four chord lengths. The minimum 
grid size in the normal direction to the wing surface is 5 x 
10“* from the leading edge to the plane of symmetry. The 
initial conditions correspond to the solution of the wing held 
at 10° angle of attack and 0° roll angle after 17,500 lime 
steps at a Mach number and Reynolds number of 0.1 and 0.4 
x I0, 6 respectively. The flowfield has no observable vortex 
breakdown. 

From the initial conditions, the wing is forced to roll to 
an initial roll angle of 8 = 10.0°. Again, by convention, a 
positive roll angle indicates that the right hand side of the 
wing when looking in the upstream direction is rolled upward. 
The wing is then released to respond to the fluid flow with a 
mass moment of inertial about the x - axis of /*, = 2.253 x 
10 r 2 Figure 27 shows the time history of the resultant motion 
and lift coefficient curve. This plot characterizes the damped 
rolling oscillations observed of wing at relatively low angles 
of attack in a subsonic freestream. At an angle of attack of 10° 
and Moo =0.1, an 80° swept delta wing will not undergo self 
sustained wing rock. From the initial displacement of 10° roll 
angle, the wing rolls to a minimum of -3.1 1® in overshoot and 
returns to a positive roll angle before reaching the steady state 
response at 8 = 0°. Meanwhile, the lift coefficient increases 
by 8.4%. 


Figure 28 shows a comparison of the Mach number con- 
tours and instantaneous streamlines of the initial conditions 
when the wing is released at 9 - 10° and the steady state re- 
sponse when the wing is at rest at 9 = 0°. Notice that there 
is very little motion of the vortex cores. As a result, the vari- 
ation of the pressure distribution is extremely small. Without 
large pressure differences between the left and right sides of 
the wing, the angular velocity remains small and the motion 
of the wing subsides. There is no noticeable lagging of the 
motion of the fluid with respect to the motion of the wing. 

Case ni-Subsonic Flow over a Delta Wing undergoing 
Self-sustained Rolling Oscillations 

The same 80° swept-back, sharp-edged delta wing of zero 
thickness is considered for this case. To duplicate the experi- 
mental investigation by Arena, (l7) the initial conditions corre- 
spond to the solution of the wing held at 30° angle of attack 
and 0° roll angle after 17,500 time steps at a Mach number 
and Reynolds number of 0.1 and 0.4 x 10, respectively. 

From the initial conditions, this wing is also forced to 
roll to an initial roll angle of 0 = 10.0° as in the previous 
case. The wing is then released to respond to the fluid flow 
with a mass moment of inertial about the x - axis of / a = 
2.253 x 1 0.' 2 Figure 29 shows the phase and time history 
of the resultant motion. From the initial displacement of 0 
s 10°, the wing oscillated in roll with a growing amplitude 
until periodicity is reached three cycles later. By i - 60, 
the motion is completely periodic with a maximum limit-cycle 
amplitude of 41.2°. For comparison, the experimental results 
for the same wing performed by Arena* 171 showed a steady 
state amplitude of 41° at the same Reynolds number Viewing 
the time histories of all three rotational properties, it is clear 
that the angular acceleration and roll angle are exactly 180° out 
of phase, while the angular velocity is nearly 90° out of phase 

Figure 30 shows the time history of the lift coefficient 
and the phase of the periodic response of the rolling moment 
coefficient. Notice that the lift coefficient curve oscillates at 
twice the frequency of the wing motion. In the phase plot 
of the rolling moment coefficient, it is interesting to note the 
three lobes of the periodic response. These lobes represent 
the energy shift from the wing to the fluid in the outer two 
lobes as indicated by the 4< +" and from the fluid to the wing in 
the middle lobe as indicated by the * 4 -”. These outer lobes 
are referred to as damping lobes. In the plot of the time 
history of the angular acceleration (Fig. 29), irregularities due 
to the damping lobes are noted near the peaks of the curve. 
Since these lobes are not present in the damped oscillation 
case, careful study of the flowfield at these points may provide 
insight into the wing rock phenomenon. 

Figure 31 shows snapshots of a complete cycle of rolling 
depicting the total pressure contours at key points labeled in 
Figs. 29 and 30. As the wing is approaching the maximum 
angular velocity, points g) to h) and j) to k), the footprint of the 
vortex core on the upward moving side appears to bow outward 
toward the leading edge of the wing. It appears that the uneven 
movement of the vortex core with respect to the leading edge 
is a result of the lagging movement of the fluid in response 
to the motion of the wing. Near the trailing edge, this effect 
is more pronounced due to the increased absolute velocity of 
the wing near the outer edges of the surface. When the fluid 
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motion catches up to the motion of the wing, the energy flows 
from the fluid to the wing promoting the rolling motion, and 
stimulating wing rock. As the wing rolls, the angular velocity 
increases until the wing exceeds 9 * ±27°. Near the trailing 
edge, the absolute velocity of the wing exceeds the limit of 
the motion that the fluid can maintain. The flowfield reflects 
this lag by the bowed appearance of the vortex core. When 
the fluid flow motion lags the wing motion, energy is absorbed 
by the fluid providing damping to the system as indicated by 
the “+" in rolling moment phase diagram of Fig. 30. As the 
wing slows, the cores appear to straighten and snap back. This 
effectively rolls the wing in the opposite direction. 

In Figure 32, a plot of the core positions at x = 0.77 is 
shown to demonstrate the symmetric motion of the vortex 
cores. Note that the cubic splines connecting the individual 
points do not represent the path taken but are merely shown 
for connectivity. During the energy transfer from the wing to 
the fluid, from points 0 to g) and i) to j), the position of the 
vortex cores exhibits a more vertical motion. When the energy 
is transferred from the fluid to the wing, the position of the 
vortex cores shift in a lateral direction paralleling the surface 
of the wing. This motion is coicident with the lagging motion 
of the fluid with respect to the wing. 

Concluding Remarks 

Recent advances in aeronautical multidisciplinary prob- 
lems have been presented. Two problems are considered in 
this paper. The first is that of fluids/structures interaction with 
applications to the twin vertical tail buffet of a generic air- 
craft. The second is that of fluids/dynamics interaction with 
applications to the natural rolling response of delta wings. 

For the buffet responses of the twin-tail configuration of the 
generic model, three sets of equations have been used for the 
aerodynamic loads, the bending and torsional deflections and 
the grid displacements due to the twin tail deflections. The 
leading-edge vortex breakdown flow has been generated using 
a 76° -swept back sharp-edged delta wing which is pitched at 
30° angle of attack. The twin tails are cantilevered without 
a trailing edge extension. Only, uncoupled bending-torsion 
response cases are considered in this study. 

The present computational results of the generic model are 
in very good agreement with the experimental data of Wash- 
bum, el. al. generic model. It is concluded that the inboard 
location of the twin tails produces the largest bending-torsion 
loads, deflections, frequencies and root bending moments when 
compared with the midspan and outboard locations. The out- 
board location produces the least of these responses. When the 
twin tails cut through the vqrtex breakdown flow, they pro- 
duce less responses due to the compensating damping effect 
produced by the left and right parts of the vortex breakdown 
flow on each tail. 

For the natural rolling response of delta wing, transonic 
and subsonic flow-regime cases are considered. The first case 
demonstrates the effects of vortex breakdown in the transonic 
regime. With the shock induced vortex breakdown, the deriva- 
tives of the motion and the aerodynamic properties show a 
very high frequency, low amplitude disturbance resulting from 
the shock-vortex interaction. Oscillations in the location of 
the shock and vortex breakdown location induces the wing to 


roll, however, the wing is unable to remain in a stable limit 
cycle. Divergence of the motion is observed when the rolling 
frequency is sufficient to cause the transverse shock to weaken 
on the upward moving side. The wing responds by continuing 
to roll until an equilibrium point is reached. 

The second and third cases are presented to provide a com- 
parison with available experimental data. Case II demonstrates 
at a relatively low angle of attack, that an 80° swept delta wing 
will not undergo self-sustained oscillations. Within one cycle, 
the wing resumes the steady state position of 0° roll angle. The 
motion of the wing and vortex cores is very slight. With the 
relatively small angular velocity of the wing, the fluid motion 
does not lag the motion of the wing. The flowfield then damp- 
ens the wing response and prevents self-sustained oscillations. 

In contrast, the third case of the delta wing at an angle 
of attack of 30° does exhibit the self-sustained limit-cycle 
rolling oscillation known as wing rock. Within three cycles of 
oscillation, the wing motion sustains a roil amplitude of 41.2° 
and a period of oscillation of 23.1 nondimensionai time. The 
phase diagram of the rolling-moment coefficient shows three 
distinct lobes which represent the energy shift from the fluid 
to the wing and vice versa. When the wing is first released 
to respond to the fluid, the pressure distribution shows a much 
stronger asymmetry than in the previous case due to the higher 
angle of attack, this results in a faster roll velocity Instead of 
the motion damping as in the first case, the increased velocity 
causes the motion to overshoot. It appears that the velocity of 
the wing near the trailing edge exceeds the ability of the fluid 
to respond. Since the motion of the vortex core is inhibited 
near the trailing edge, the core exhibits a distinctive bowed 
shape. While the motion of the fluid lags the wing response, 
energy is stored in the vortex cores. The wing motion slows 
as a result of the damping provided by the energy transfer to 
the fluid and reverses roll direction. With the slowing of the 
roll rate, the motion of the fluid ceases to lag the motion of 
the wing. The vortex cores appear to snap back. On doing so, 
the energy stored in the fluid is imparted to the wing causing 
the roll velocity to increase. This cyclic motion builds until 
the energy transfer of the system is balanced and the periodic 
response is maintained. 
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Figure I. Three-dimensional grid topology of the twin tail-delta wing configuration (Midspan) 



Figure 2. Three-dimensional and front views showing the total pressure on the surfaces, and 
the vortex-core streamlines. Initial conditions (Inboard position). 



Figure 3. Initial conditions for static pressure and instantaneous streamlines in a cross-flow 
piane. x = 1.096 (Inboard position). 
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Figure 7. History of the deflection and load responses for an uncoupled bending-torsion case, 
co = 0.3, a = 30 °,R e = 1.25xl0 6 , (Inboard position). 




(b) Right tail 


Figure 8. Total structural deflections and root bending moment for an uncoupled bending- 
torsion case .V/oo = 0.3 , a = 30°, R t = 1.25x10 s , (Inboard position). 



Figure 9. Three-dimensional and front views showing the total pressure on the surfaces, and 
the vortex-core streamlines. Initial conditions (Midspan position). 













Figure 10. Initial conditions for static pressure and instantaneous streamlines in a cross-flow 
plane, x = 1.096 (Midspan position). 
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Three-dimensional and front views showing the total pressure on the surfaces, and 
the vortex-core streamlines. Uncoupled case after it = 9,600 (Midspan position). 
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Figure 12. Snap shots of static pressure and instantaneous streamlines in a cross-flow plane, x 
= 1.096. Uncoupled case after it = 9,600 (Midspan position). 




(a) Left tail 


(b) Right tail 


Figure 13. Distribution of the deflection and load responses for an uncoupled bending-torsion 
case. Moo = 0.3, a = 30 °,R C = l.‘25xl0 6 , (Midspan position). 


(a) Left tail 


(b) Right tail 


Figure 14. History of the deflection and load responses for an uncoupled bending-torsion case. 
Moo = 0.3, a = 30°,f? e = 1.25xl0 6 , (Midspan position). 




(a) Left tail 


(b) Right tail 


Figure 15. Total structural deflections and root bending moment for an uncoupled bending- 
torsion case. M, 0 = 0.3, a = 30 °,R e = 1.25xl0 6 , (Midspan position). 




















(a) Left tail (b) Right tail 

Figure 18. Distribution of the deflection and load responses for an uncoupled bending-torsion 


case. .V/ 0 c = 0.3, a = 30 °,R e — 1.25x10®, (Outboard position). 
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Figure 19. History of the deflection and load responses for an uncoupled bending-torsion case. 
Moo = 0.3, a = 30 °,R e = 1.25xl0 6 , (Outboard position). 




Figure 20. Total structural deflections and root bending moment for an uncoupled bending- 
torsion case. Moo = 0.3, a = 30°, R e = l.‘25.rl0 6 , (Outboard position). 


History of Lift Coefficient History of Drag Coefficient 




Figure 21. History of the forces on the wing/twin-tail configuration. Moo = 0.3, a = 30°. R t = 
1.25zl0 6 , (Inboard position). 
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History of Lift Coefficient 



n 


History of Drag Coefficient 



Figure 22. History of the forces on the wing/twin-tail 
1.25xl0 6 , (Midspan position). 

History of Lift Coefficient 



configuration. A /,*, = 0.3, o = 30". R f = 


History of Drag Coefficient 



Figure 23. History of the forces on the wing/twin-tail configuration. M 
1.25x10 , (Outboard position). 


= 0.3, Q = 30°. R, = 


Transonic Flow-Initial Conditions 
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Figure 24. Mach Contours near the Wing Surface with Surfaces of Constant Entropy (s 
= 0.5) and Instantaneous Streamlines; a = 20.0°, 0 = 0.0°. * 
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Transonic Flow-Divergent Oscillations 




Tune History of Roll An* * -Angular = 0.0^ - 
and Lift Coefficient; Moo = 0.&5 ,R t - ^ jIlu ’ a 
n It (noints of interest annotated). 




Figure 26. 


Snapshots of the Mach Contours hear the Wing Surface with Surfaces of Constant 
Entropy (s = 0.5) at Points of Interest at a)0 = -0-47 . b) 

0.44°, d ) 9 = 11.08°, e) 6 = 24.27°. 
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Subsonic Flow-Damped Oscillations 
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Figure 27. Time History of Roll Angle, 6 , Angular Velocity, 9, Angular Acceleration, 0 and 
Lift Coefficient; M x = 0.1, R e = 0.4x10 6 ,q = 20.0 °.0, c = 10.0°. 0, r = 0.0 “//. 


9 = 10 . 0 ° 


9 = 0 . 0 ° 
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Figure 28. Comparison of Mach number Contours and Instantaneous Streamlines of Initial 
Conditions ( 0 = 10.0°) and Steady State Response (0 = 0.0°). 

Subsonic Flow-Self-Sustained Limit-Cycle Oscillations 
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Figure 29. Time History of Roll Angle, 0, Angular Velocity, 9, Angular Acceleration, 9 
and Phase of Angular Velocity, 6\ A/^ = 0.1 . R t = 0.4j10 6 .o = 30.0°. 0, r = 
10.0°, 9, c = 0.0 °jt (with points of interest annotated). 
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Figure 30. Time History of Lift Coefficfient and Phase of Rolling Moment Coefficient indicat- 
ing Energy Transfer Lobes, (with points of interest annotated). 
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Figure 3 . . Snapshots of a Complete Cycie of RoHing Osciltaion depicting te Tonsure 
Contours at Points of Interest at 0 9 = 41.1 . g) ^ - 27.3 , h) 0 0.0 . 

-40.8°. j) 9 = -27.5°. k) 0 = -0-2°. 0 0 = 41.2°. 
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Simulation of tail buffet is studied for several delta wing-vertical tail configurations. 
Flow conditions are chosen such that the wing primary- vortex cores experience vortex 
breakdown and the resulting turbulent wake flow impinges on the vertical tail. The 
dimensions and material properties of the vertical tails are chosen such that the deflec- 
tions are large enough to insure interaction with the flow, and the natural frequencies 
are high enough to facilitate a practiced computational solution. This multi-disciplinary 
problem is solved sequentiedly for the fluid flow, the elastic deformations and the grid 
displacements. The flow is simulated by time accurately solving the laminar, unsteady, 
compressible, Navier-Stokes equations using em implicit, upwind, flux-difference splitting, 
finite volume scheme. The elastic vibrations of the tail are modeled by coupled bending 
and torsion beam equations. These equations are solved accurately in time using the 
Galerkin method and a five-stage, Runge-Kutta-Verner scheme. The grid for the fluid 
dynamics calculations is continuously deformed using interpolation functions to smoothly 
disperse the displacements throughout the computational domain. Tail buffet problems 
are solved for single tail cases, twin F/A-18 tail cases and twin highly swept generic 
tail cases. The use of an apex flap for buffet control is also computationally studied. 
The results demonstrate the effects of inertial structural coupling, Reynolds number, aft 
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fuselage geometry and spanwise tail location on the tail buffet loads and response. Favor- 
able comparisons with experimental data Indicate that the present aeroelastic method 
is well suited to providing qualitative insight into the tail buffet problem, as well a$ 
quantitative data for refined long duration simulations. 



